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In the present work, experimental and analytical investigatioris 
have Been carried out on some dynamic prohlems in beam -foundation 
interaction. Effect of size, shape of the footing, degree of f lexihil.' ■ 
of the footing and magnitude of force levels to which the footing is 
subjected, have been experimentally investigated. Various types of 
soil media, such as a loosely compacted, uniform, medium to fine dry 
sand} fine sand stabilized by adding sone percentage of cement* and 
a natural silty soil have been selected for the studies. 'Prom the test 
results, it has been found that the maximum displacenent varies 
proportionately with the force applied. The effect of the flexibility 
on the spring constant was found to be negligible as predicted from 
theoretical results. On the effect of the size, and shape of the 
footing, the spring constant ratio and dacgjing factor showed variation 
with length/width ratio. The spring constant ratio have been found 
to be increasing with D/b ratio while the damping factor is decrerstr; 



with L/b ratio. Jin atten^t has Iseen made to analyse this observed 
behaviour. 

Solutions for the response of beans carrying concentrated mass, 
resting on soil medium have been presented idealizing the soil medium 
as Winkler model and as Pasternak model. I)an5>ing has also been included 
in these models. The beam with concentrated mass attached to it and 
resting on elastic foundation has been analysed for steady state dynamic 
force. The displacement magnitudes have been calculated for infinite 
and finite beams with different values of concentrated mass attached 

i 

to it. It has been found, from the analysis, that with concentrated 
mass increasing, the maximum displacement at the centre and at points 
along the length of the beam, decreases. The rate of decrease in 
maximum displacement have been found to be nonlinear with concentrated 
mass. The displacement magnitudes have also been evaluated at various 
points along the length of the beam, for beams of different length. Prom 
the theoretical results, it has been observed that along the length 
of beam, the maximum displacement decreases, with the rate being non- 
uniform along the length of beam. However, with the beam length, the 
maximum displacements have been found to be increasing, though marginally. 

Beams resting on elastic foundations subjected to impulse loads 
have also been analysed. A general solution for displacement response 
at any point along the length of the beam has been derived through 
numerical inversion of Laplace transforms. Post andWidder's formnla 
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has been used for numerical inversion of Laplace transform for infinite 
beam on elastic foundation, Berger's technique of Laplace inversion 
for infinite and finite beams has been used to get the numerical values 
of the response functions. From the results, it has been observed that 
with damping, the maximum displacement decreases, ilso, along the 
length of the beam, the maximum displacement has been fo\ind to be 
decreasing. 
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CHAHER 1 


HffiBOHiaillOF 


1*1 GtEETEBAIij 

Ihe dy naml o a n alysis of foundations finds many apiSLications in 
modem tecshnology sudi as strip footings, railway trades, aircraft 
XTinways, buried pipes, rocket testing tracks etc, Sich problems of 
foundation-soil interaction are generally solved by incorporating 
the reaction from the foundation into the response mechanism of the 
structxire by idealising the foundation by a suitable mathesDoatieal 
model. In majority of the cases, the msponse of the structure at 
the contact surface is of prime interest* Ihus it would be of 
immense help in the analysis, if the foundation can be represented by 
a simple mathematical model with the parameters of the model diaracte3>- 
izing the true behaviour of the system as closely as possible, ■ To 
accomplish this objective, many foundation models, from the M.hkler 
foundation to the elastic continuum idealisation have bee® evolved. 

A comprehensive review of some of these foundation models is given in 
reference (29) . 

for apy method of analysis, the critical step would be 'ttie 
correct evaluation of s<ril parameters -ghich characterize the 
bdiaviour of soil m^a. These soil parameters wiJi be different 
for different media* Rirthermore these pajameters will depeai on 
loading conditions, st^rain distributions developed in the soil mass^ 
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size and shape of the st 3 cucture etc, Ihroug^i various testing 
procediires, either by testing the soil samples in the laboratoiy or 
by performing the tests in~situ , these paiamete3?s are evaluated. 

In both types of testing, the values of these parameters depend on 
testing procedure, environmental conditions and the extent of 
medium, !Ehen the question arises: Can these values be used as 

such for the analysis of various foundation soil interaction probleims 
in all situations? 

!Qie dynamic loading effects on stationary structure have been 
the subject of great interest, A lot of worir has been done on the 
bdiaviour of beans subjected to steady state forces and transient 
loads. But for beams, resting on soil media, the available literatTxre 
is quite meagre. In technical applications it is not infret^ent that 
a beam resting on soil media idiidh is under vibration carries a lumped 
or concentrated mass somewhere along the length of the beam, A 
machine resting on a deep beam is an exampLe, Ho solutions have been 
reported so far on aich jproblems also, 

!Che vibrations, due to aircraft landing on airstrips, due to 
rocket movement on rocket-sled-tradc, due to tremsieiri: loads on buried 
pipes, due to dynamic loads on colxunn-foundation system etc. are few 
examples which require Ihe analysis of beam-foundation i^stem, simulating . 
the field conditions, !Qie dynamic soil paredneters, also, regfutlre to 
be critically examined Aether th^ reqp;dre any ioodiflcations in their 
values for drajsg^d situations of the vibrating systems. Some theoretical 
and experimental investigations are pr^ented and analysed for scmae 
aspects of these- poroblems. 
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1.2 BBIEP BByiEff OF LCTERAJOBE 

The earliest fonnolation of elastic foundations v?as due to 
THiirikler vdio assumed that the foundation model consisted of closely- 
spaced independent linear springs (29). The itthkler model approxi- 
mates the reaction of the continuum by a simple expression vdiich 
results in a single ordinary differential equation governing -ttie 
responses of the beam. The solution is strei^t-forvirard for any 
loading configuration. A more difficult analysis by Idealising 
idle foundation as an elastic continuum or models derived from the 
conbioium hypothesis is ofdy possible in most cases through essentially 
approximate and mmerical methods. Several other ^thois have evolved 
various partly continuous models conceptually de-veloped from the WLrKLer 
model| such as I^stemalCy iilonenk;o~SOrodi<±y W.egiardty. Heteryiy 
Kerr etc, (29). The objecti-ves in employing a model are mathematical 
simplicity and relati-vely accurate representation of the real gystem, 
Ofteuy the Visco-elastic model ( Spring-Bashpot system) is used in 
dynamic analysis ovilng to math eRiiit tool simplicity^ thougi the above 
mentioned models ihich are used static loads, can also be used for 
dynamic analysis with appropriate inclusion of damping. In visco- 
elastic model, Idle soil medium is Idealised in terms of piysical 
constant s-suoh as spring constant, damping coefficient and the equi-valent 
mass of soil participating in the -vibration. Attempts have been made 
to evaluate, these parameters appropriate to represent the media as 
closPly as possible. Tletcher and Hermann (l7) dev^oped foundation 


■# 
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peirameters based on elastic constants, consideilng the foundation 
model in iidiich the soil reaction is given by, 


2 4 

/ . 3 w , 3 w 


( 1 . 2 . 1 ) 


#iere q(x) is foundation reeiction, w the deflection, and t the 
foundact ion modulus. In another discrete model (54), the foundation 

3 ?eaction has been expressed in tenns of spring constant and damjing 
constant as follows 


q(x) = (k^+ ia^c^) 


( 1 . 2 . 2 ) 


inhere q(x) 


m 


- foundation reaction 
= static stiffness of the footing 
s= dynamic spring factor 
ss dynamic dampLag factor 

r 

UJ • 

= nondimensional frequency, _ 

s 

= circular fredpiency 


r = radius of the equivalent circular tooting 
eq 

V = daear wave velocity in the medium. 

3 

These spring and damping parameters and c^ are the functions 
of coefficients that have been calculated and presented for various 
fourdatlon chapes by Veletsos and VejSjic (55) for footings resting on 
elastic half apace, on 1iie basis of Jhimelee^B (38) work. &ico and 
Westmann (3l)» Bainer (41 ) Sarraain, Boceet and fliitman (47) pre^nbed 


the variation of the dynamic spring factor pasjaaeter with Ihe fre^enc^ 
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of vibrations. Stiucture-Jfoundation infeeraotion studies were extended 
by Jennings and Bielak (24), Novak (35), Saulad^nski (52) etc, for 
various other aspects of loading conditions on the basis of this soil 
reaction* Analysis of the change in the dynamic spring factors and 
the damping factors were carried by Eosset, Whitman and Dobjy (45) 
and ifeiiner (42). Ehe damping is attributed to radiation of energy 
due to the propagation of the waves in the medium, and the energy 
dissipation throu^ intergranular friction. Ihrou^ experimental 
evidence, Kausel and Bosset (27) indicated that the energy dissipation 
takes place largely due to internal friction. It was also indicated 
that this energy dissipation is independent of frequency. About the 
spring constant, most of the methods presume that the spring constant 
is independent of amplitude of motion, OSiis assumption is laot strictly 
valid (59). However, so long as the motions are within certain limits, 
the behaviour of foundations msy be taken as linear. 

Prom the experimental and theoretical work, it is dear that the 
stress distribution for the static load at -aie interface of structxuje 
and soil may be quite different from the presumed distribution in the 
Wirkler model. Yet from a bending moment diagram of an infidte beam 
acted upon by a concentrated load, derived by Biot (s), it is evident 
that for certain values of spring constant based on Mriler model, these 
bending moment diagrams are in reasonable agreement (56). After 
conducting some experdmeiis, it was shown by Vesic (57) that WirKLer 
hypothecs is pjractically satisfiad, not only for lor^ bease but dso, 
for beams for moderate lengths, ^ the beams of infinite length on 
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an ideally elastic confcimum are arjalysed by wi rik-T er foundation, 
the bending moments are overestimated, contact pressxn:e and 
deflections of the beam are underestimated (56). 


Using Biot’s theory of bending of beams, Vesic (56) derived 
the expressions for the coefficient of subgrade laaction K ^ , K , 
for an infinite and finite beam respectively developed a 
relatioru^p between them in terms of elastic constants as given by 


K = t.23 


B . B 
s 


LOpCl-u^) E^lJ 1H. 


0,11 


B 


s 




K 


Tidiere 



modulus of elasticity of Beam 
I =» moment of inertia of Beam 
Eg= modulus of elasticity of soil 
UgS® mission’s ratio of soil medium 
G^sa coefficient of pressure distribution 


( 1 • 2, 3 s) 


(l.2,3b) 


(1.2.4) 


Vesic and Johnson (58) substantiated later the reliability of the ei^iresedons 
by model studies of beams on silt, A varisction of subgrade reaction 
with i/B ratios is given in Bowles (il). 
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Prom the theory of beams on an elastio foundation a foundation 
can be evaluated as ri^d or flexible. Heteryi (20) proposed the 
flexibility criteria depending on the value of XI (vvhich considers 
the width, length and elastic properties of the nediam) , where 



in udiich K is modulus of subgrade reaction, 

3 


Vesic ( 56 ) also proposed the criteria that was quite similar 
to the one proposed by Heteryi. Depending on the flexibility of 
the beam, Vesic reconnnended the procedure for approximate analysis 
(lable 1.1). - 

If the foundation is absolutely flexible, ^loaded by a vertical 
pressure, then the stresses in the soil under the fouMation will be 
distributed uniforEily, but settlement under the foundation will vaiy. 
!I!he average settlement is directly proportional to the pressure applied, 
with the coeffiaLent of proportionality ' O^, called the coefficient 
of elastic uniform compression, given by 




(1.2.6) 


Using Sciileidier* s solution (5) for an average settlement, a 
coefficient 0^ vhich depends on l/B ratio, can be derived. Then 
* C^* can be written as 


G 


u 




■I 


.(1.2.7) 
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where is coefficienb of elastic uniform compression, 
Pg is uniform pressure under the foundation, 
is average settlement under the foundation, 

A|j is base area of foundation 

E 

0 ^ is ^ v,E elastic constants 

r ^ d. ^ 

1-v 

She expression for C is given by 

s 

ir VT 



(1.2,8) 


where a a l/B 

Ihus the spring constant K for absolutely flexible fourdation is 
given by 

O,*’® (1.2.9) 

(l-») 

there Ii is the length of the foundation. 

Itor an absolutely rigid rectangular foundation Gorbunov-Posadov (3) 
computed the values of the spring constants and it was found that 
there is hardly a difference of 3 ^ in the values, A comparative values 
of Cg for different i/B ratios for rigid and flexible foundation 
are presented in the table 1 . 2 , A general reviet? on the ^namic 
soil-structure interaction is provided by Eameewra Eao (26), 
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theoretical investigatioias on beam vibrations find its origin 
in curvature-bending reiationMiip derived by Bernoulli, Beiiiaps 
Euler and Bernoulli were the first to find the fundamental frequency 
of a cantilever bar and apparantly, Ludwig was the first to deal with 
steady state dynamic solution of an infinite beam on an elastic 
foundation (32). In 1954 Kenny (28) included linear viscous damping 
in the basic equilibriina equation and solved for the corresponding 
stea<ty state response, the steady state solution with a forcing 
function cos mt instead of 1 was obtained by Mathews (33)* 
Hbwacki (36) solved the basic equilibrium equation. 


4 2 

EI-“+Kw+m a E 6(x). 6(t) (1.2.IO) 

3x 3t 

Eor the taransienfc loads, Sfcadler and 2ireeves (50), Hades (40) 
presented the analytical solutions for infinite and finite beams 
a:espectively for the dyi^amic response of the Euler-Bemoulli beam 
including the effects of damping, axial load or shear layer etc, as 
given by the equation 



oil + 
° at + 


= q(2:>t) 


(1.2.11) 


■vdaere I is axial load or shear coefficient 


Based on Galerfeen's method, Eringen (16) deadLved deflection 
equations for transverse, time dependent, force, Boley and Chao 
(9,10) analysed, approximately the TimosheEfco beam (Botary inertia, aM 
daearing force of the beam is included in the basic equation), based 
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on. a ^Izavelling Wave* approach and the principle of virtual woik, 
M L k lowitcz (34) in. his analysis of fTimosiherfco beam, divided the 
total deflection into its bending and daear components and solved 
these eqvdLvalent set of coupled equations. Anderson (1) presented 
the general series solution for the flexural vibration of Timoshenko 
beam. lor dhort beams, some discrepency in the deflection from 
theoretical ^ear coefficient of a Timoshenko beam was reported by 
Kri^namurthy (30) after conducting scane experimental investigations. 
Thus he jxresented a new fomnulation in vhich governing equations are 
two simultaneous partial-integro-differential equations. 

lor beams on elastic foundation, some analytical results are 
available vhere the foundation modulii are nonlinear (l 9 » 18 ). 

Salto and Tomizo (46) studied the effect of mass of the soil in deriving 
frequency equations for an infinite beam resting on elastic foundation 
using the wave propagation theory, Berger ( 5 ) considered the effect 
of ^ear force, rotary inertia and viscoelastic damping on the 
response functions of Ihe beam, 

IQae problem of lateral vibration of a simply supported beam 
with concentrated laass attached at the centre was studied by 
Hoppman II ( 22 ) , taking the Euler-BernouUi beam. Chen (12) 
pxreserted a new fomulation in which the dirao-delta function is used 
in differential equation, Irequeaacy equation for the no r mal modes 
of vibration were obtained by Baker (2), Scimth and Ifes ( 49 ) 
pn^esehted the analysis for simpQy supported beam having rotary inertia., 
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later Das and Deidxmuldi (l5) derived a general expression for 
Timosheifco beam with concentrated mass at the centre. So far, no 
solution is reported for beams on elastic foundations with concen- 
trated mass attached to it. 

On the effect of impact on the beam resting on the elastic 
foundation, there is very little information available. Besides 
ITowacScL (56) , Hoij)iaan (21) derived expressions for simily supported 
beams with or without damptiiJg, using energy method. Schwieger (48) 
collected some experimental data to check the central deflection of 
the beam, Prasad (39) also conducted some experiments a>id model 
analysis for theoretical prediction of response for impact load, 

1*3 SOOPB OP THE PBESERT IWESlIGAflOH 

In this thesis experimental investigations have been carried 
out to study the effect of size, shape of the foundation, degree of 
flexibility and magnitxide of force levels to iidiich it is subjected. 
Attempts have been made to analyse these observed results in the 
li^t of available theories on dynamic soil behaviour, 

She present woifc also includes the study of beams 37esting on 
soil media suh;jected to steady state vibrations and impulse loads, 
Ihe soil media is idealized as WirScler model as well as jfestemak 
model with damping included in both the models, A beam with conceit, 
trated mass attached to it, is analysed for steady state vibrations, 

A general solution for displacement at ary point along Ihe lei:^& of 
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■fctie beam subjected to impulse load is deiived throu^ numeidcal 
inversion of laplaoe transform for finite and infinite beams on 
elastic foundations. 

ffltie procedure results of experimenbs to evaluate the soil 
parameters throu^ cyclic plate load test, cross-Ji ole— seismic testj 
block test are presented in Chapter 2. Erperiments have been conducted 
to study 1iie effect of size, ^ape of the structure, degree of 
flexibility of the structures and magnitude of force levels on these 
soil parameters. Chapter 2 also presents the experimental details 
arws results of these studies. On the basis of existing theories on the 
b^aviour of soil media an attempt is made to analyse these results. 

To veri^ the Winkler model hypothesis, a flexible beam resting 
on ideal springs is tested arid the 3 resulte are compared wiiaa the 
analytical solution for sinusoidal force. 

In chajter 3, the beam carrying a concentrated mass and resting 
on soil medium is analysed for steady state vibration. Two mathew 
matical models, namely the Wiriaer model and the ^stemak model with 
and without demping are taken for the analysis of aiLer'-iBernoulli 
beams* Scpressions arid results are presented for the general solution 
of finite and infinite beam s with concentrated mass at the centre of 
the beam. These solutions TfUslch c»n be used, for any boundary corsiLtions 
are checked wi1±i the exLstijrg sohitlojH for the beams carrying a 
concentrated mass with end supports (without foundation). Variation of 
magnitude of dlsiiLace3saents along the length of tiie beam for various 



13 


values of concentrated mass, shear values of Iksternak foundation 
model, and with the change in the length of the beam is presented 
in this chapter. Comparisons have been made with some existing 
solutions. iixperimental results for the response of an beam resting 
on soil medium with concentrated mass attached to it have also been 
discussed* 

For £inite and infinite beams, resting on kLastic foundation, 
general solution for displacement at any point along the length of the 
beam for an impulse load have been derived Ihrou^ numerical inversion 
of the lapLaoe transform of the basic equations of motion. In 
caaapter 4, various inversion tedaniques of Laplace transform are 
enumerated and Best & YELdders^t and Berger' s techniques are used in 
deriving the solution. Solutions have been compared with existing 
results, 

Ihe general conclusions and the scope for further work have been 


presented in Chapter 5 
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1.4 TO Allots 

Each symbol ia defined when first introduced and some of the common 
symbols used in this thesis are listed below 


A 


A = 

c 

A 

max 


^1 

B = 


b 

C 

a 

d 

a 


shear coefficient of Pasternak foundation 
Base area of foundation 
cross sectional area of the beam 
maximum amplitvjde 
non-dimensional frequency 

s 

A^,A^,Ag = arbitrary constants 
width of the beam 



shear parameter 

laplace transform of ^ear parameter 


C 

r 


E 


1-v 



a coefficient 

coefficierri; of elastic uniform ccmipression 

initial constants 

damping constant 

critical damping constant 

dynamic damping fhctor 

external daaplng 

internal damping 

coefficient of pressure distribution 



D 

a 

5 

a 

E 

E 

s 

e 

o 

E 

F 

o 

g 
I 
j 


= damping factor 
ss shear parameter 

= lapiace transfoim of ^ear parameter 
= modiilus of Elasticity 
•1» elastic modulus of beam 
= modulus of elasticity of soil 
= coefficient of restitution 
= a concentrated sinusoidal forcing function 
= force of constant magnitude 
= shear modulus of soil 
s= arbitrary constant 
= moment oi|* inertia of beam 
,= Bessel's function 


- 1 / 4 ^ 1 / 4 " 

= spring constant of the soil medium 
= coefficient of subgrade reaction for finite beam 
= coefficient of subgrade reaction for infinite beam 
=s modulus of sub grade reaction 
= static stiffness of footing 
= tynamic spring factor 
k^jkgjEj^ foundation model parameter 
L = lengtii of the foundation 
iCe) =s lapl^e transform of the function P 

M = concentrated mass attached to the beam 


K 

K 

E 

K 

K 


st 


M(z) = bending moment cxf beam at z 
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m = unit mass of the beam 

= mass of falling body 

= arbitrary constants 
n s= integer number 

(Ph»P2) 

« Jacobi Pciynomial 

P(z) = soil reaction on beam at z 

(p) = transfoim variable 

p = uniform pressure under the footing 

s 

p^jPg = arbitrary constants 

q(x) = foundation reaction in vertical direction 

k 

q(x,t) = foundation reaction in vertical direction at a distance x 
at tine t 

E = nondime nsional beam length = l/ X 

as 

^ ~ as 

n 

r^^ = radius of the equivalent circular footing 

= Average settlement under the footing 
S(z) = Slope of the deflection curve at z 
=: oxLal load 

1D_ s= time of contact 

ij 

t(z) s= shear force of beam at z 

t = time parameter 

V = shear wave velocity in the medium 

s 

V = velocity of freely falling mass just before collision 
o 

W = wei^t of falling body 

w as deflection in vertical dire^ion 



= foundation deflection beyond the length of the beam 
= distance paiameter 
= non-dimensional vertical displacement 


y 

y^ 

Kp) 

y(2) 

y(z,0) 


z 

a 

6 

e 

Y 

o 

y 

y 

0) 

(I) 

n 

e 

♦ 


W 

= — = non-dimensional displacement of foundation b^ond beam 
X 

= laplace transform for y(z) 

s= non-dimensional vertical displacement at z = w/X 

= i'yirw^jinfln.cf i nnft.1 vertical displacement at a distance z and 
at time 0 , 

= non-dimensional distance a/X 

= concentiated mass ratio = y and l/B in equation ( 1.2.8) 

c 

= i^y 

“ S' 

r= 4(l+2I>p + P^) 

_ 1/2 

- V (j ; 

= arbitrary constant 
= mass density of soil nEdia 

= 4Al-r^)^+ (2Dr)^ 

= Ibission’s ratio of soil medium 
= circular frequency 
=5 mtural circular frequerry 
=s nbrwddmensional time variable 
= jhase angle 





X 


X 

fi(t) 

n 


= Caudiy’s function 


= Characteristic length = 



^^2 .4^4xl/2vl/2 

= (2m + 2(m +n ) ) 

=Kroneckex delta of x 

ssKronecker delta of t 

= characteristic shear -value of soil medium 


Ip s= a function 

5 = notation . 
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Table 1,1 

Criteria for Hexibility of Beams 


S,lIo, 

Size of Beam 

Criteria for 

Procedure 



Distinction. 

Recommended 

1 

long beam 

XIi > 5.0 

Winlsaer analysis 
assuming infinite 
beam 

2 

Moderately long 

2,25 < Xl* < 5.0 

WinlfLer analysis 

3 

Saorb moderate 
beam 

0.8 < XI < 2.25 

De Beer Method 

4 

Saort beam 

'll < 0.8 

Treat as perfectly 
rigid 


(Based on Yesic 56) 


Table 1,2 

Ysilues of Cg for flexible and rigid beams for various 

I/B ratios 


l/B 

1 

1.5 

2 

5 

5 

10 

16 

C for flexible 
beam 

1.06 

1.07 

1.09 

1.13 

1.22 

1.41 

1.507 

C for rigid 
® beam 

1.08 

- 

1.10 

1.15 

1.24 ’ 

1.41 

mm 


Barfcan ( 3 ) 




CHAITEE 2 


EXlEEIMimAIi IHVESTIGATIOUS 


2« 1 GEHBEAIi 

Por the expeadjaercfcal iiwestigations on dynamLc soH^stnict^ire 
interaction, a variety of soil media; a loosely compacted, well 
graded, medium to fine sand; a fine sand-stabilized by adding some 
percentage of cement; and a natural soil deposit are taken, iUiese 
investigations are carried out to study the effect of daape, size of 
footing on the foundation parameters of soil, represented by a 
Visco-elastic model (spring and dash pot in series, that is Winkler* s 
model with viscous damper)* Studies also include , the effect of 
magnitude of force level to ifdaich the footing-soil system is subjected 
and the effect of flexibility of footing. Experimental results of 
a beam resting on ideal springs are compared with analytical 
solution of a beam resting on a Winkler foundation. In this chapter, 
the procedures and results are presented for these investi^tions 
besides the soil properties and its evalxiation procedures* 

2*2* 1* Soil Broperties-Byaluatioa Procedures, 

!I!he procedure for basic ihysioal soil properties are quite 
standard and are available in jnost of liie text books on sdl 
mechanics, Por dynamic studies, soil parameters in teaans of ^ear 
modulus, Poisson* s ratio, elastic modiaus or alternatively, spring 
comtant and damping factor of the soil medium need to be evaltrated. 
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IS 5249-1969 ( 23) provides a standard pa?ocedure for obtaining these 
parameters, 3y conducting cyclic plate load test, the soil 
parameters for dynamic studies can also be ascertained in terms of 
coefficient of subgrade reaction and coefficient of elastic urdfoim 
compression (3) while from the block test (23), the soil parameters 
are evaluated in terms of elastic constants. Another procedure, 
Cross-hole test, was presented by Stokoe II and Woods (51) which 
facilitates the evaluation of the dynamic soil parameters in terms 
of elastic constants. Ihis procedure seems to be better, compared 
to block test, as the time taken in setting up the experiment and 
conducting the test is far less, and the ability of this method to prob' 
deeper depths. In evaluating the parameters, however, all the 
three methods were adopted. TShile the block test procedure is 
quite established (25) , the cross-hole test and cyclic plate load 
test are described veay briefly, 

2«2,1,1, Cross-hole Method. (Stokoe II & Wood) 

Ihe procedure followed is to advance two vertical boreholes 
some distance apart into the ground to the same depth, A transducer 
is placed at the bottom of one of the bor^oles, She transducer 
is a vertical, velocity pickup attadaed to a pipe that allows the 
transdvuser to be placed at the bottcm of the borehole and can be 
removed. An impulse rod is placed down tt^ other boa^ole, A 
small top plate with a vertical velocity transducer attadaed at the 
bottcm of the first borehole is connected to the input of an 
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oscilloscope while transducer attached to the impulse rod is connected 
to the trigger of the oscillascope. 


Ihe test is conducted by striking the plate on the top of 
impulse rod with a hammer which triggers the oscilloscope and also 
sends a compression wave down the rod, GDhe compression wave from 
the rod starts the wave propagation into the soil at the bottom of 
the rod. Body craves are generated in the soil and their arrivals 
at the velocity transducer in the opposite borehole are disiieyed 
on the oscilloscope, This procedure of striking the impulse rod 
and recording of the wave arrival is done several times with the 
equipment in the initial position to determine an average time, illhe 
impulse rod and velocity transducer are removed and placed down 1±ie 
opposite borehole and the test repeated. Then the seismic 
equipment removed from the bor^ole and boreholes are advanced to the 
next depth and the procedure is repeated. The pxrirary wave 
velocity and diear wave velocity are determined by dividing the 
distance between the bordhoies by total travel time taken by primary 


wave and shear wave respectively. Using the density valu^, the 
shear modulus, elastic modulus and Itoisson's ratio are evaluated. 



can be written as 


Baissou’s ratio 
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Tdiere v is primaiy wave velocity in the medium 
0 

v„ = shear wave velocity in the medium 
s 

p =mass density of the medium, 

2,2, 1.2 Qyclic Hate load Test (Barkan 1962). 

In this method, a square or circular rigid plate is loaded 

and unloaded in incremercbs, with the record of . settlement and 

recoveny in settlement, in the soil beneath the plate, IDae jlate 

is reloaded to a hi^er increment and then urCLoaded, Bach hi^er 

loading increment makes the plate settle more and unloading allows 

the plate to reduce this increased settlement* Ihe part of the 

settlement Tdaich does not recover after unloading is called 

permanent settlement and the remaining part of the settlement which 

is recoverable on unloading, is called elastic settlement, Ihis 

elastic settlement inereases with increase in vertical stress, Hie 

ratio of the vertical stress to this settlement is called the 

coefficient of elastic uniform compression. In the present work, a 

2 

circular plate of area 1640 cm was used for the test, Slg, 2. 1 
shows the load and total settlement (elastic + plastic) relation. 
Hois relation is quite linear, Results of cyclic plate load test 
are ^own in lig. 2.2, Out of the elastic and plastic settlement, 
the elastic settlement is separated out and plotted with vertical 
stress in lig, 2.3, It can be observed from the graph, that the 
elastic settlement is proportional to pressure applied. 
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2# 2# 2 Soil Eroperfcies**Results , 

Por the investigations, a variety of soil media were takeru 

In the laboratory, a well graded dry sand, with a uniformity coeffi»* 

cient of 2»1, was filled in a large container# This saM with 

3 

density of 1*62 gey cm and void ratio of 0«65 was poured 
compacted in layers of 10 cm throu^out the area by a rammer 
imparting the equal energy everyvfaere# In the field, two types of 
soil media were selected* Ifihile one soil medium was a natural 
soil deposit of silty sand (Density 1,62 gn/cm?), having liquid 
Lhnit aM plasticity index equal to 27,5 aM 10 respectively, the 
other was micaceous fine sand filled in a trendi of size 20m x 5m x 2in, 
Ihis sand was stabilized by adding 5^ cement and compacting it at 
the optimiim moistiire conbenb, !llie demity of the soil was 1,55 gir/cm^. 

^toile the soil in the container and in the trendi was more 
.homogeneous, isotropic and uniform, the natxiral soil deposit has 
a range of particle sizes, !Ihese three different soil media were 
selected to represent various soil conditions, She loose diy sand 
may absorb the wave ener®r more thus increasing the damping coefficient, 
the stabilized soil may increase the stiffness of the soil thus 
increasing the spjring constant, Hfhile both these soils are extending 
to a l i m i t ed length, the natural soil deposit extends upto infinite 
extent. 

Evaluated soil paroperfeies are presented in Stebles 2,1, 2.2, 

2# 3* Ihere Tables 2, 1 and 2,2 fibow Idie p£Eima:Qr, di^r wave velocity, 
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elastic constanfcs, soil demity etc, for stabilized and. natural soil 
deposits respectively, Table 2,3 presents the summary of these 
elastic constants and soil parameters for all the soil media, 

2.3 IgHAMIO SOIL-STHJCTUIiE DCTEMCTION STUDIES 


The effect of size, shape of the footirag, magnitude of force ■ 
applied and the flexibility of the footing on soil -structure intera- 
ction are dealt in this section. Experiments were conducted on 

(i) a circular plate of 22.9 cm diameter and thickness of 1,25 cm, 

(ii) a square plate of size 20.4 cm x 20,4 csm x 1.25 cm, (iii) a 

rectangular plate of size 40.8 cm x 10.16 cm x 1,25 cm, and (iv) a 

long strip of dimension 81.3 cm x 5,08cm x1,25 cm, resting on three 

different soil media. The base area of all the plates and beam have 

2 

been chosen to be same and equal to 413.0 cm , with the wei^t approxi- 
mately 4.0 kg. The frequency independent sinusoidal force was applied 
on the centres of these footings by an electromagnetic exciter. The 
magnitude of the force level was var i ed on these footings, EoUowicg 
sub-sections describe, the experimental set up, instrumentation, 
sequence of operations, results and discussion of the results. 

2.3.1 Experlmertal Set-up. 

The various instrumerdis used for excitation, measurei»nt ani 
control of dynamic characteristics of Ihe expjeilmental system are 
briefly presented. Gare was taken at esuh step to check the cali- 
bration and rating of the prcimary equipcffints during their operation. 
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Eada reading was repeated several times in the course of single 
measurement or over a period of time to ensure the consistency of 
the results. 

2,3. 1.1 Test Erame for Ijynamic Jests for Soil in the Iiaboratoiy. 

A schematic diagram of the fraiae is shown in 5ig, 2,4. It 
consists of a large bin filled wi-th loosely compacted, well 
graded sand, !I3ie plates and beam were tested one after another, 
after keeping them on the surface of sand so that unifoim contact 
between soil and the footing is maintained. Before, during and 
after the testing, the footings were checked for tiniform contact of 
these with soil, and horizontal level of these plates, The exciter’s 
moving table was fixed at the centre of the footings. To keep the 
exciter in inverted position so that its moving table ^ust jrests on 
the footings without adding any weigat of its own on the footing, 
a frame made up of steel angles was fabricated, Bie position of 
the exciter could be vertically adjusted by means of a central 
screw arrangement in the steel frame. The frame is big enougi to 
accommodate the sand container. The exciter was fix^d ti^tly to 
the frame so that at the time of resonance, thei*e i^ould not be any 
mDvemenrt: of the exciter. The moving table of the exciter was along 
the axis of the footing. 

2*3* 1.2 Test Rrame for Byxiamic Tests of Soils in the Sield , 

Slg, 2,5 presents a schemetic diagrame of the frase, for testing 
in the field. This frame made up of steel angles, was fabricated 
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with modifications suitable to field conditions. Instead of a 

cenbral screw arrangemenfc , the exciter in a steel cabinet was 

kept in inverted position by steel ropes. The wei^t of the 

exciter and the cabinet was balanced by three sets of equal 

wei^ts, hanging on the same rope on the other end as ^own in 

the sdhemetic diagram (lig. 2.5), The rope passes throu^ a set 

) 

of frictionlesa pulleys specially fabricated for the purpose. 

During preliminary investigations, yhen the exciter was kept by a 
central screw arrangement, it was found at certain frequency, 
this frame was vibrating at its resonant frequency. Thus, by 
adopting this method for keeping the exciter in inverted pjosition 
a n d lowerirjg it to the footing so that it just rests on the 
footing, the problem of frame resonance was avoided. Besides Hie 
fact that the exciter does not add its own wei^t to the footing, 
another advantage derived from this method, was that during the 
vibration, footing did not loose its contact with soil medium at 
ary time, 

2.5»2 Inst rument at ion 

2,5.2, 1 Oscillator 

Jin oscillator (ifodel 200JIB Audio Oscillator) capable of gemra.- 
ting sinusoidal pulse of constant voltage in the fret^enty range of 
20-.10000 ops was used for excitation. Maximum outprt power of the 
signal was 1 watt with an overall accuracy of 2^, The desired 
frequency and output voltage were controlled by the frequency dial 
and amplitude knob xespjectively. The output signal was amplified 
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2,3. 2. 2 Amiaifier. 

A power amplifier (Model 2120 MB Power Amplifier) used to 
amplify the out jut signal from the oscillator, Ihe frequency of 
sinusoidal pulse remained unaltered and only the output TOltage was 
amplified. Hie signal distortion throu^ the amplifier was less 
than ^ within the frequency range of 20-1000 cps, The maidmum 
output power of the amplifier was 125 TA and input signal voltage 
was 1 V, mns, 

2.3.2.3 fedter . 

The exciter used was an eiectiromagnetic device (ifodei Efl. 1250 MB) 
■srfaich could produce the maximum force levels of 12 kg and maximim 
displacement of 1,25 cm at different pgy loads and accelerations. 

The exciter had a frequency range of 5-10,000 cps with total wei^t 
of moving table 0,523 kg. Its signal distortion was less than 5^, 

Hie frequency of the exciter was controlled throu^ the oscillator. 

The maximum power input to the exciter was limited to 15 volts, 

7 amps, 76 watts and controlled throu^ an amplifier. 

2.3. 2.4 Yibmtion Transducer, ( OB-1566 P-52). 

In the accelerometer type of piotaip which was used in the 
exi>6iimenfcs, the transducer is mde of a piezo-electric crystal, 

A piezo-elect lie crystal developes an electrical charge directly 
proportional to the pressure applied to it, Hiis pressure is due 
to the acceleration of the base on which the transducer is fixed. 

This ji<dEup was used to measure amplitudes of dispiao^eut, velocity 
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aM acceleration of Idae base, Ihis was held in position with the 
magnetic clamp or fiaed otherv/ise. 

2. 3. 2.5 Yibiation Meter (ge-1553 AK) 

iEhis instrument consisted of a sensitive vacuum toibe voltmeter, 
iEhe vacuum tube voltmeter included a hi^ gain band amplifier, a 
caliberated attenuator, average reading, peak reading, and peak to 
peak reading voltmeters, integrating circxiita to convert acceleration 
signals to velocity or displacement signals. Its operation range was 
to 2 - 20,000 cps. An outlet was provided to connect it with the 
oscilloscope. 

2.3 .2. 6 Oscilloscope. 

itor the visual display of waveform, a cathode ray oscilloscope 
with storage facility was used for visual cheek of forcing function 
or alternatively to facilitate to photograph the waveform trace. It 
consisted of a cathode rey tube in which a beam of electror® Ulvuni- 
nates a spot on the screen, Ifhis spot can be accurately positioned 
both horizontally and vertically by application of voltage across the 
pair of deflection plates between which the beam passes. One set is 
used to control the vertical position and the other,horizord:al position 
of the beam, Por most applications, a calibrated saw-tooth voltage 
is applied to the horizortal deflection plates so that beam moves from 
left to ri^t at a known veiocHy and rapidly returns to the left side 
of the tube to repeat its sweep from left to ri^t. 2iis provides 
a time scale on the screen of calhode ray tube, Kis ^gnal to be 
measured is amplified and applied to the vertical deflection plates. 
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2.3,3 Sequence of Operation. 

Befoie starting the test, the soil, contained in the bin, 
and in the field, on viiich these jlates or beam rests, was made 
even and level. It was ensured that footing is in contact with 
soil medium at every point and was level, Ihe axis of the moving 
table and the axis of centre of gravity of footings were checked 
to coincide and remain vertical, The level of the footing and 
the axis of vibration were checked intermittently. 

After resting the footing (circular, square, rectangular 
plate or beam having the sasae base area of 413 cm ) on the soil 
media, the exciter was fixed on the footing. The exciter was 
connected to an oscillator throu^ an amplifier. An electric 
sinusoidal force signal was fed from the oscillator to the power 
amplifier from viiere, after getting the signal amplified, it was 
fed into the exciter. The exciter, in turn, converted this 
electric signal into mechanical oscillation vdiich was trarumaitted 
to the footing soil system. As mentioned earlier, the frequency 
and the amplitude of force level was controlled throu^ the oscilla- 
tor. The nature of simsoidal force produced by exciter was 
visually verified thixju^ an accelercanfiter attadied to exciter, which 
fed the signal to oscilloscope throu^ a vibrationmeter. 

She response of the footing soil system was measured thrc«j#i 
a system of vibration pickup and vibrationmeter. The displacements, 
velocities and accelerations were measured at various points on the 
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footing throu^ an accelerometer type of pickup TAfaida was attached 
to the footing, Ihis transducer picked up the vibaration and fed 
into the vibrationmeter, Ihe nature of oscillations of the footing 
soil system was monitored on oscilloscope by taking out the signal 
from -vibrationmeter and feeding it into oscilloscope, iig, 2.6(a) 
shows a block diagram of the test set up. Starting from -very low 
frequency, -the responses were noted in tems of frequency, dispiace- 
ment, -velocity and acceleration of the footing soil system over a 
large frequency rarge, OIhe force level was constant throu^out 
the frequency range aM was increased for the next set of readings, 

For each type of footing (circtilar, square, rectangular shapes or 
beam) three to fi-ve force levels were applied for the response curves , 
ae frequency ranged from 300 to 600 cps. 

2,3«4 Bresenfcation of Besults . 

Ihe experimental results, of -the tests conducted on circular 
plate, square plate, rectangular plate and long strip (beam) resting 
on three different soil media, ara presented in -this section, lEho 
base areas, the thickness and the wei^ts of these footings were 
the same, Sable 2,4 provides relevant infoimation of these footings. 

2.3«4,1 lests for &?il in the Laboratory, 

In -the first set of experiments, -these four footings (circular 

plate of dia 22.9 caa, square plate of 20.4 can x 20.4 cm, rectangular 

plate of 40.6 cm x io.2 cm and a long strip of 81,2 cm x 5.08 cm)of 

2 

base area of 413 on , resting on the soil contained in bin were tested 
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i3i the laboratoiy^ !Ehese were subjected to sinusoidal f 03 ?ce at 
their cent3?es. The magnitude of force level ranged from 3^0 kg 
to 9*0 kg# The response oulves were obtained for deflection;, 
velocity and acceleration at each frequency. The velocity, acce*- 
leration values wei^e utilized for the verification of dispLacemenfc 
values* 

Ei-g* 2#7 ^ows the response cuarves of a circular plate for 
four different force levels# The plot is betv/een frequexacy (cps) 
versus magnitude of displacement (in m#m#)# It can be observed 
that resonance occurs at 385 cps in the frequency range of 300--600 cps# 
With increase in force level, the resonant frequency does not cdaange# 
Only inaxunum deflection is increasing with inciease in force level# 

OSae frequency dispELacemenfc curve for sqxiare plate is shown 
in ELg# 2#8# me response curves for four different force levels 
were obtained as for circular plate# 

For this plate also, frequency ranges from 500 to 600 cps with 
the resonance at the frequency of 385 cps# With increase in force 
level ,maxinixim displacement was increasing v^iile resonance frequeiiKy 
was not mch differejit# 

For rectangular plate, the resonance shifted to a hi^er 
frequenxy 410 cps* ELg. 2*9 ^ows 'itxe response curv® for four 
different force levels* liagnitude of displac^ent increase with 
the increase in force level# For "fee rectangifLar plate, the 
displacements we3?e measured at more points than for circifliar plate 
or square plate aM an average value was plotted with the freqaen<y# 
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In 5ig^ 2.10, frequeccy-displacemertt response is pioftted for 
the long strip. During testing, it was observed that the dispLace- 
ment along the length of the beam was varying at any particular 
frequency. Ihis figure shows the deflection of the beam at its 
centre. It can be seen that like other footings, the resonant 
frequency was not varying much with the magnitude of force level ard. 
remain at 5 10 cps. Only TtiayfrinTn deflection is increasing with force 
level. 

Prom these figures (lig. 2.7, 2,8, 2,9, 2, 10), it can be 
observed that resonant frequency is increasing from square plate 
to beam. As the base area, thickness, and the mass of the footings 
are same, the increase in fre(jiency and Ttia.-iri'muin displacement, could 
be attributed to the ^bape of the footing, Wa-Hmum displacement 
and frequency are plotted against l/B (the length of the plate or 
bear/widlii of the plate or beam), in Jigs. 2.19 ard Sig, 2,20 for 
all the soil media. Another observation can be made frcaa Pigs. 2,7, 
2,8, 2,9,' 2, 10, that there is some proportionality in the increase of 
mariiium displacement with increase in force levels. Both the marl. mum 

displacement and resonant frequency are increasing with the l/B 
ratio as diown in (Pi^, 2,19, 2,20), 

Prom these curves, the spring constant and the damping factor 
are calculated assuming the plates and the beam as rigid, lable 2.5 
presents these values for various force levels applied on them with 
the formulae used in Ihe calculatiois to evaluate these values. 
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2,3»4.2 lest Eesults for Sfcabilized Soil Medlijm* 

The response curves for circifLar plate, square plate, rectangular 
plate and long beam resting on stabilized soil are presertbed in 
Slgures 2.11 to 2.14 respectively. Stor this soil medium also, 
fair to five force .levels, varying from 3.0 kg to 7.0 kg vjere selected 
for all the footings, Ihe frequency ranges from 300 cps to 600 cps. 
like in previous case, the displacements for rectangular plate were 
measured at more points and average was plotted (ELg. 2,13) vhile 
for long beam, the dis placements at the centre of gravity were measured 
several times and plotted in ilg, 2, 14. 

iOiere was hardly ary variation in resonant frequency with 
magnitude of force level. With i/B ratio, the maxiiiiim displacement 
and resonant frequency were increasing (ELg, 2.19, 2, 20), as was 
found for the soil in the container. Same w^, there was some 
proportionality noticed in the increase of maximum displacement with 
the force level. 

treating the plates and beam as rigid, the spring constant 
ard damping factor are calculated and tabulated in Table 2,6, 

2,3 .4.3 Test Results for Matrtral Soil Deposit. 

In the last set of results, the circular plate, square plate, 
rectangular plate and the long beam, resting on the natural soil 
deposit, were subjected to stea<3y state oscillation of constant 
force as was the case for the previous two typ©3 of soil media. 

2, 15 to 2. 18 '.abow the frequency' displaceraent curves for three 
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to four force levels. As usual, the frequeocy ranges from 300 ops 
to 600 cps. Maximum displacement increases with force levels as 
can be seen from the iigures (2.15, 2,16, 2,17, 2,18), Maximum 
dispLacemeiit also inc3?eases vaLth i/B ratio, Ihe resonant frequency 
does not vaiy much with foarce levels. It varies with l/B ratio 
only (iig» 2*19, 2.20), Table 2,7 presents the spring constant and 
damping factor values for these footing. These values are 
calculated treating the plates and the beam as rigid. 

Taring the averages of the values of spring constant and damping 
factor, calciJlated from the response curves for circular, square, 
rectangular plates and the long beam on the three different soil 
media, the variation of spring constant and damping factor with i/B 
ratio is plotted. Pig, 2,22 shov^ the spring constant ratio 
variation (spring constant for ai'^r i/B ratio/spring constant of 
circtiLar plate of the same base area) with h/B ratio. The spring 
constant ratio is increasing with increase in l/B ratio. Pig, 2,23 
ptresents Idae variation of damping factor with increase in l/B ratio. 

The damping factor decreases with the increase in l/B ratio. The 
rate of increment of spring constant ratio, aai rate of decrement of 
damping factor is not the same for three different soil media, hut 
■Haere is a definite ti«nd with increase in l/B ratio, 

2*3*5 Hscussion. 

Based on these observations, the effect of magnitude of force 
levels, the effect of size acd shape of the footing, and the flexibility 
of the footing on dynamic foundation parametea® are discussed in this 


section* 
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2.3.5. 1 Mggaitude of Force. 

Irom the response curves (Hgs. 2.7 to 2.21 and Table 2.5 to 2.7 
it can be seen that there is some proportionality in the increase of 
maxLmna displacement with the increase in force levels. Generally, 
the increase in displacement is linear with force levels, as in the 
theoretical analysis the maximum displacement is directly proportional 
to the force level, ahus, it may be taken that the maxiaxam displace- 
ment is directly proportional to the force level. 

2.3.5. 2 Effect of Size and Siape. 

2.3. 5. 2.1 Spring Constant . 

In the theoretical visco-elastic model, the footing, resting 
on it, experiences same stiffness at all poijts below the footing 
viien the footing-soil system undergoes vibr^ation. In practice, 
however, it is not always true. Stresses below the plate vary from point ^o 
point depending on the nature of soil medium and flexibility or 
rigidity of the footing, A rigid or flexible plate on soil does 
not have uniform stress distribution. It may vary from csiter edge to 
Itaier area, Also with the change in the base area, if the load 
distribution remains the same, the stress below the footing at a pcaidi 
dees not remain constant (3). Thus liie coefficient of elastic uniform 
compression depend on Ifae area and the shape of the foafeng. In 
another words, the spring constant will depend on the shape of tSae 
ooting, if the base area remains the same. 
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Below a circular or square plate, the stress distribution 
and thus the stiffness at the edges will be different from that 
in the inner portion. Ihus the net sum of the stiffness wil3. be 
less than the values if tho stiffness had been iinifoim throu^out 
the area below the plate. Having the same base area, if liie length 
inci^eases and width decreases, the stress distribution becomes more 
and more uniforn, IDhis uniformity in stiffness brings the net 
sum to a little hi^er value. (Hius, as the length increases and 
width decreases, the value of spring constant appears to be increasing, 

Por ri^d rectangular plate, Gorbunov-Posadov calculated the 
values of a coefficient vixich can be equated to spring constant ratio 
for various values of l/B ratio, Barkan ( 3 ) has presented those 
values in tabular form. According to them, if the base area of a 
footing remains constant, the coefficient of elastic uniform compre- 
ssion increases with an increase in i/B ratio. Ihe variation of 
spring constant ratio based on this theory is shown by dotted line in 
figure 2.22. 

2.3* 5 .2. 2 Damping Pact or. 

Prom the figure 2.23, it can be observed that damping factor is 
decreasing with the increase in l/B ratio, axis damping liiich the 
soil medium possesses, arises because lixe ener^ is dissipated throu^ 
intergranular friction and due t o wave pxropagation in the medixm ( 42 ) • 
She energy dissipation throu^ intergrranular friction is called 
material damping and reveals in a hysteretic load deformation curve 
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for the soil. She energy loss due to propagation of waves is called 
geometiie damping or radiation dampirg, Shis damping is a function 
of frequency of the system (59) while the material damping or internal 
damping is frequency independent to a lai’ge estent (27). 

G3ae decrease in damping m^ be attributed to the increase in 
stiffness and probably the reduction in material damping. As the 
damping factor is a function of damping constant, spring constard; and 
mass of the soil participating in vibration, the increase in spring 
constant m;^ result in a decrease in the damping fact or. As analysed, 
it is tafcen that with increasing l/3 ratio, the stiffness of the soil 
will increase, thus, resulting in the hi^er spring constant value. 

2. 5. 5*3 Effect of ^legibility. 

Based on elastic theory, Barkan (3) has presented the relationship 
between spring constant ratio and l/B ratio (Equation 1.2.8) for 
flexible rectangular plate and compared these values with the values 
for rigid plates* She difference between these values is not larger 
than 3 percent. Prom the experiments 1iie value of spring oonstant 
ratio is calculated for the long beam with reference to sqcP-re plate, 
treating the beam and square plate as rigid in one case, and as absolutely 
flexible in the other case, lable 2.8 presents the comparison of the 
these two cases, Bor this i/B ratio, the theoretical values of spring 
constant ratio for flexible and rigid cases are also presented in this 
table, it<mi the experimental results, the spring constant ratio does 
not seem to be different for flexible case in comparison to rigid case, 
■ffiithin experimental errors, these values be taken as e<|uel to the 
theoretical values. 
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2*4 ICTSBIMEHrAJIOISr lOE WIICLBE MOIiai VBEIMGAIIOI?’ 

To sfcudy the response of beam on elastic foundation, often the 
foundation is idealised as ISCLricler inodel. An experiment has been 
designed to veii:^ the theoretical values based on TiELEkler model, 

A schematic diagram is shown in liguro 2*6, 

2.4.1 Experimental Setup 

Beam of mild steel and of unifoim cross-section was put on 
ten almost identical linear springs at a spacing of 10, l6 cm, !Btie 
details of the characteristics of the beam and springs used for the 
experimental study are given in table 2. 9, These springs were welded 
to a fairly heavy daannel making the base rigid. The ends of the 
beam were kept free, Thus a uniform free-free Eule 2 >-Bemoulli beam 
resting on IKlrkler model of foundation was simulated in the experi- 
mental set up. 

In studying experimentally, the response of the beam subjected 
to steady state sinusoidal force, various parameters such as displace- 
ment, acceleration, velocity of the beam vrere measured, Instrumeriation 
arrangements were the same as for the studies on dynamic footing-soil 
interaction. Block diagram is ^cwn in iig. 2,6(a), IiHtiumentation 
details are given in the section 2,3,2, 

All the ten spring which were used in the model were tested on a 
uni— axial compression testing machine. It was found that the ic^ 
deflection relationship w^ fairly linear during loading, tihioading 
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and reloading. Table 2.9 presents the value of spring constant, 
diaract eristic length etc. 

2.4»2 Results and Discussion . 

Ror three different set of force levels, the beam was subjected 
to steady state sinusoidal foi-ce vtlch varies from 5.15 kg to 12.0 kg 
Prequency ranges from -30 cps to 120 cps. The displacement, accele- 
ration, and velocity were measured at centre of the beam and along the 
length of the beam at the interval of 10*2 cm, A typical set of 
measured displacements along these points are presented in Table 2,10. 
y/ith the increase in frequency, the maadjmjm dispLacemenb was increaising, 
and it reached its first peak at the frequency of 55 cp>s, indicating 
that beam is under first mode of vibration. In this mode, mfl-riTmi-m 
displacement occurs at the centre and the point of zero deflection, at 
l/4th distance from the one end of the beam. Besides the certre of 
the beam, displacements were hi.^er at the free ends, 

Erequency displacement responses at the centre of the beam, for 
the force levels 5,15 kg, 9,4 kg and 12.0 kg are shown in figure 2,24, 
leak of the displacement at first mode of vibration occurs at 55 cp^. 
The variations of displacement alor^ the length of the beam for few 
values of the frequencies of vibration are presented in Hgures 2.25, 
2,26, 2.27 for force levels of 5.15 kg, 9.4 kg, and 12 kg respect! veiy* 
Irom these figures, it can be observed that the amplitude of deflection 
is decreasing at eveiy point on the beam as the frequency goes on 
iaacr^slng. Also, the pcxint of zero deflection diifts towards the free- 


free end 
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Based on the theoretical solution of a beam resting on the 
Winkler foundation, the values of the amplitude of displacement 
are calculated. The expression is given by 

(l+cos pR cosh hr) (cos ]LZ + co^ ]jz) 

(sin yR cosh yR + cos UR sirh uR) 


y(z) = - 


R X 


2EI ^^3 


(l+sin UR sihh yz) (cos yz - cosh yz) 

+ ' ■ ■ “ ■" ' ■ 

(sin uH cosh UR + cos uR sirh uR) 


(2.4.1) 


Ror the first mode of -vibration of beam, the amplitudes of 
displacement along the length of the beam, for all the -three force 
le-vels were calculated and are plotted by dotted lines in 
Jigures 2.25, 2.26 and 2.27. 

The experimental curve is in good agreement with the values, 
calculated from theoreticaQ. analysis of Beam on Wirfcler foundation. 
This probably reinforces the confidence in modelling the foundation 
medium as per M.rfcler*s hypothesis. 
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lable 2. 1 

Eropeities of Sbabilxzed Soil Medium 


Bore ’ 
log 

Depth 
in m. 

Time 
taken 
by P 
wave in 
millisec. 

Time 
taken 
by S 
wave 
in 

millisec, 

P Wave 
velo-. 
city 
n/sec. 

• 

S Wave 
velo- 
city 
n/sec. 

Poiss- Den- 
ion' s sity 
Eatio gin/ 
cc. 

E 

kg/ 

cm^ 

Siear 

Modulus 

kg/cm^ 

G*Om 
Stabili- 
zed soil 

0 

39 

64 

256 

156 

.30 

1.60 

1030 

397 

sand 

1.83m 

0.60 

39 

63 

256 

159 

.30 

1.55 

1040 

400 

Silt 

clay 

2.13m 

Silt 

clay 

1.22 

38 

63 

263 

159 

.35 

1.55 

1090 

400 

with 

Karikar 

1.83 

38 

62 

263 

I6l 

.35 

1.55 

1090 

410 


2 2 
Distance between Trigger hole Average G=402 kg/cm Average E=1060 ’kg/cii' 

and Eeceiver iiole=lOia, Poisson's Eatio=0,32 Desnity=i,56 gn/c, c. 
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Table 2.2 

Ero parties of natural Soil Deposit 


Bore 

Depth Time 

Time 

P wave 

S wave 

Dois- 

Density 

E 

& 

log 

in m taken taken 

velo- 

veloc- 

son’s 

gn/cn^ 

Elastic 

Shear 


by P 

by 

city 

ity in 

Batio 


modtilus 

modulus 


wave 

S wave 

w/ sec 

n/sec. 



kg/cm? 

kg/om^ 


in 

in mil- 








milli 

li sec- 








sec 









0 . 00 “ 








Silt 

with 

0 33.0 

51 

280 

180 

0,250 

1.60 

1275 

509 

clay 

0.62 31.5 

50 

291 

183 

0.266 

1.65 

1425 

562 


2.13m 








Silty 

clay 

1 . 2 C 30.0 

50 

305 

183 

0.380 

1.60 

1520 

548 

with 

Eatasnr 

1.83 28.0 

48 

327 

191 

• 

1.65 

1790 

615 

3.36m 








Sandy 

silt 

2.44 29.0 

49 

315 

186 

0.370 

1.65 

1620 

585 


mtii 

Bajari 

“ ' ' ^ ■ ■ ■ ■■ ■■■' ■ ■ ■ ■■■■ II IT ■ -- 11^ 

Average G=550 Kg/cjii Possion's ifeitio=0,31 
Distance between trigger bole and 2 3 

Eeceiver bole=9 , 15in, Average B=1460 Kg/cm Density = 1.62 @i/cai 
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Table 2.3 

Suamaiy of I^roperties of Soil Media on vdiich the Experiments 
Tsrere conducted 


Type of Soil Medium 

Medium and line 
Uniform Sand in 
the bin 

Stabilized Sand 

natural Soil 
Deposit 

Soil identification 

Badari«r Eine 
uniform sand 

Micaceous sand 
stabilized with 

5^ cement 

Sil-ty clay 

Void Eatio 

0.60 

0.65 

0.50 

liquid Limit 


- 

27.5 

Elasticity index 

- 

- 

10.0 

Total unit wei^t gn/cc. 

1.62 

1.56 

1.62 

Total mass density 

, 2 
kg sec 

cm:^ 

1,65 X 10~^ 

-6 

1.59 X 10 

1.65 X 10"^ 

2 

Elastic modxiLus fc^cm 

- 

1060 

1460 

2 

Shear modulus kg/cm 

395 

402 

550 

BaLssion’s Batio 

0.30 

0.32 

0.31 
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lable 2.4 

Data regarding footings used in the Experiments 


Serial Eroperties 
Eo. 

Circular 

Hate 

Squ&re 

Hate 

Rectangular 

Hate 

long Strip 

1 DLidb nsions 

22.9 cm dia 
(9'0 

20.4 x2Q.4cm^ 

(8" X 8”) 

> 

40,6cm X 
10.2cm 
(l6’'x 4") 

81,2cm X 5.08em 
(32" X 2") 

2 (Thickness 

1.27 cm. 

(0.5 '0 

1.27cm. 

(0.5”) 

1,27 cm. 
(0.5") 

1,27 cm. 

(0.5") 

3 Wei^t 

4.048kg 

4.055kg 

4.032kg 

3.93kg 

4 E 

ModTilus of ELa^ 
sticity 

6 

2.1 X 10 
kg/cm2 

6/2 6 
2.1 X 10 kg' cm 2.1x10 

kg^cm^ 

6 , 2 
2.1 X 10 kg' cm 

5 I 

Moment of Iner- 
tia 

- 

3.48 cm^ 

4 

1,74 cm 

4 

0,87 cm 
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Table 2«5 

Data for Soil Contained in bin, emd foimilae used for 

evaluation 


a, QircxiLar Hate 


Elate dimension 9” dia.l/2” thickness Wei^t = 4*048 kg. 


Sr.Hb. 

Force in 
leg 

Besonant 
frequency 
in cps. 

Maximum 
Amplitude 
in sum. 

Spring constaid: Damping factor 
in kg/cm. 

1 

3.0 

380 

0.0072 

24000 

.090 

2 

5.1 

385 

0.0100 

246CX) 

.104 

3 

7.0 

385 

0.0110 

24900 

.129 

4 

7.5 

385 

0.0220 

24500 

.072 

b. Square Elate 





Elate dimension 8" x 8" 

X 1/2" 

wei^t = 

4.055 kg. 

Sr.IIb. 

Force in 
kg 

Besonant 
ftequenqy 
in cps. 

i^ximiua 
amiiitude 
in Bum. 

Spring constant lamping factor 
in kg/cm. 

1 

3.2 

385 

0.0075 

24600 

.090 

2 

5.0 

385 

0.0090 

24800 

.113 

3 

5.3 

390 

0.0095 

25400 

.110 

4 

7.0 

385 

0.0180 

24550 

.080 

c, Bectangia.ar Elate 




Elate dimension l6"x 4” 

X 1/2” 

Wei#^it = 

4.032 kg. 

Sr.IIb. 

Force 
in kg. 

Besonent 
fre(gfienicy 
in cps. 

TnaTTfirtiim 

ampiitiide 
in Bum. 

Spring conartant* Damping factor* 
inkg^cm. 

1 

3.1 

410 

0.0055 

27909 

.102 

2 

5.0 

410 

0.0099 

27800 

.100 

3 

5.2 

415 

0.0180 

28000 

.051 

4 

7.5 

420 

0.0225 

28400 

.058 


* Assuming the' jOLate is rigid 
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Sable 2.5 (confcd.) 
d, long Beam 

Beam Hmensioa 32" x 2” x l/2'* Weigjit =5.95 ki 


Sr.Ko. 

Porce 
in kg. 

Eesonant 
frequency 
in ops. 

Islaximum 
amplitude 
in m,m. 

*Spring constant 
in kg/cm 

Damping factor''* 

1 

3.4 

500 

.0133 

39620 

.032 

2 

5.0 

510 

.0174 

41230 

.035 

3 

7.4 

510 

.0230 

41260 

.039 

4 

9.5 

520 

.0270 

42910 

.041 


*Asstiiiiing tHie plate as rigid 


e, Pormulae Used 


iferameters 


BOnnulae 


Known parameters - P, f . max, m,P = Porce, flft=Eesorient frequencgr = 


fn = natxizal frequency 


fa /l-.2ir . A = Maximum Amplitude.msmass of Plate 
7 max 

1 

2ir m- 


I) = Bampiag factor 


1-21) /- ii^.ms *^>2 


2D05r^ 


-) (2 fffto)' 


K = Sparing constant 


■Snax. 2D 1-D^ 
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Table 2.6 

Data for Stabilized Soil 

a. Circular Hate 

Hate dimension = 9" dia l/2" thickness wei^t of the plate 

4.048 kg 


Sr. Ho. 

Itorce 
in kg. 

Resonant 
ib?equency 
in ops. 

Maximum 
Amplitude 
in m.m. 

Spring constant Damping factor 
in kg/ cm. 

1 

3,0 

375 

0,0062 

23900 

.105 

2 

4.2 

375 

0.0075 

24200 

.123 

3 

5.3 

380 

0,0086 

24600 

.129 

4 

6.0 

385 

0.0100 

24900 

.122 

b. 

Square Hate 

Hate dimension = 8" 

* 

X 8" X 1/2” 


wei^t of the plate 
= 4.055 kg. 

1 

3.3 

385 

.0075 

24500 

.090 

2 

4.6 

385 

.0078 

25000 

.120 

3 

6.2 

385 

.0096 

25100 

.130 

4 

7,1 

385 

.0100 

25100 

.130 
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Table 2.6 (confcd.) 


o, Eectangia.ar Hate 

Hate dimension = 16*' x 4” x l/2” Y/ei^t of the ia.ates4.032kg. 


ar.M-. 

Porce 
in kg. 

Besonanb 
Prequency 
in cps. 

Maximum 
amilitude 
in m.m. 

Spring constant 
in kg/cHu* 

*Bamping pactor 

1 

3.1 

390 

.0070 

25000 

.089 

2 

5.6 

400 

.0115 

26200 

.060 

3 

5.0 

400 

o 

LTV 

O 

• 

26250 

.064 

4 

5.1 

395 

.0110 

26000 

.091 

5 

7.2 

395 

.0125 

26100 

.110 




^Assuning the beam as 

rigLd 

d. 

Long Beam 






Beam dimension = 32” x 

2” X 1/2” 

Wei^t of the Beam 3.93kg, 

Sr.N. 

Porce 
in kg. 

Besonant 
Pre^enc^ 
in cps. 

Maximxim 
amplitude 
in m.m. 

Spring constant *DampiDg Pactor 
kg/cm.* 

1 

5.1 

490 

0.0130 

38000 

.010 

2 

4.5 

495 

0.0185 

38700 

.010 

3 

6,1 

500 

0.0223 

39500 

.011 

4 

7.5 

500 

0.0290 

39500 

.011 


^Assuming the beam as rigid 
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Sable 2.7 

Data for Hatural Soil Deposit 


a. Circular BLate 

Slate dimension = 9" dia. l/2” thickness Wei^t of the ia.ate=4.0481i:g. 


Sr.lTo. 

Force 
in kg. 

Resonant 
Frequency 
in cpe. 

Mazinium 
amplitude 
in Huzn. 

Spring constant 
in kg/cm. 

Damping factor 

1 

3.0 

390 

0.0066 

25000 

.090 

2 

4.0 

390 

0.0079 

25050 

.100 

3 

5.5 

395 

0.0096 

25800 

.121 

4 

6.2 

390 

0.0103 

25150 

.119 


Square Hate 





Hate dimension 8" z 8" 

X 1/2" 

Wei^tof the 

plat e=4. 055kg, 

SuHo. 

Force 
in kg. 

Resonant: 
frequency 
in CPS. 

Mazimum 
amplitude 
in m.m. 

Spring constant 
in kg/cm. 

Damping factor 

1 

3.1 

385 

0.0070 

24500 

.090 

2 

3.9 

390 

0.0079 

25400 

.099 

3 

5.0 

390 

0.0096 

25500 

.103 

4 

6.5 

395 

0,0116 

26000 

.108 



c* Eectaogular Elate 

Elate dlmenaion = i6» x 4” x l/2» Weight of the ia.ate=4.032 


St^O» 

Force 
in kg. 

Besonant 

Frequency 

in CPS, 

Maximm 
amjQ.itude 
in rcum. 

Spring constarcb 
in kg/cm,* 

Lamping factor* 

1 

3.3 

445 

.012 

32250 

.043 

2 

5.5 

425 

o 

m 

29700 

.060 

3 

6,7 

425 

.016 

29800 

.065 




*Assuming the beam as ligid 

d. 

Long Beam 






Beam dimension, = 32" x 

2" X 1/2" 

Weight of the 

ja.ate=3.93 

Sr.Hb, 

Force 
in kg. 

Besonant 
Frequencgr 
in cpe. 

MaTriTmiTTi 

amplitude 
in m.m. 

Spring constant 
in k^cm,* 

Lamping factor* 

1 

3.2 

485 

.015 

37300 

.0090 

2 

4.3 

490 

.018 

38600 

• 

o 

o 

3 

5.1 

490 

.022 

38500 

.0100 

4 

6.3 

495 

.026 

39200 

.0099 


♦AsstiiiiiiJg the beain 6.s rigid 


1.1. I. ¥ r4^rsM 

Aoc. Ml. A 
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Table 2.8 

Increase in Spring Co ret ant Ratio for l/B=l6 


i/B = 16 

Treating beam Tareating Beam as Beam as 

as rigid beam as Rigid ELexible 

flexible 

1.683 1.810 1.50 1.51 


1.595 1.505 1.50 1.51 

1,520 1.432 1.50 1.51 


Table 2.9 

Details of Beam and Spring used for HrfeCLer Model Idealisation 


Beam Sa.ze Densilgr Young' s 2 Spirii:^ 

length Width Thickness in kg^cn? modulus Kg. cm L constant K 

in can. in cm, in cm, of ESLas- kg/cm? 

ticity 

kfi/cm2 

101,5 5.04 1.26 7.8x10”^ 2.1x10 1.76x10 3.2 7,0 


Beam on Sandy Soil 

Beam on Stabilized 
Soil 

Beam on latural Soil 


Experimental Values 


Theoretical Values 
Barkan (3) 
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Table 2, 10 


HLsilacemerit Values with Srequency at R>rce level 9.4kg. 


- — 12 !' 
— 8 " 



16 



A. 



Fo siA/ cat 

B F 


Sr.lfo. 

Erequ- 

eaqjr 

cps 

Disila- 
cement 
at A 
in m.m. 

Ihspla.- 
cemenfc 
at B 
in m.!n. 

Us pla- 
cement 
at C 
in num. 

Displace- 
ment at 

D in m,m. 

DLs place- 
ment at 

E in Bum. 

DLsplacemeat 
at E. 
in num. 

1 

35 

0.800 

0.680 

0.540 

0.400 

0.300 

0,050 

2 

40 

0.900 

0.720 

0.440 

0.200 

0.032 

0.200 

3 

45 

0.750 

0.550 

0.260 

0,020 

0. 170 

0.200 

4 

50 

0.800 

0.520 

0.220 

0.130 

0.350 

0.400 

5 

55 

0.700 

0.460 

0.100 

0.240 

0.500 

0.550 

6 

60 

0.450 

0,280 

0.024 

0.220 

0.350 

0,400 

7 

65 

0.400 

0,220 

0.028 

0.190 , 

0.300 

0.350 

8 

70 

0.220 

0.120 

0.035 

0.160 

0.230 

0.200 

9 

75 

0.160 

0.075 

0.046 

0.140 

0.190 

0. 180 

10 

80 

0.150 

0.050 

0.050 

0.110 

0.150 

0.150 

11 

90 

0.060 

0.022 

0.027 

0.064 

o,oeo 

0. 100 

12 

100 

O.C^ 

0.027 

0.037 

0.076 

0.080 

- 

13 

110 

0.066 

0,022 

0.030 

0.063 

0.067 

- 

14 

120 

0,063 

0.015 

0,022 

0,050 

0.045 







Pressure (kg/cm 












Steel 

frame 


Vibration 

pickup Natural or stabilized soil 


Fig. 2.5 A schematic diagram of the test frame 
for field tests 
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Fig. 2. 6 (a) A block diagram showing the arrangements 
of instruments for constant amplitude 
excitation on the footings 



E Sin U)t 
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To vibration meter 


Ideatized Foundation Model 


Fig. 2v6 (b) tetealized Winkler Model 








Frequency (Cps) 

Fig. 2. 7 Displacement frequency curve for soil contained 
in the binCfor circular plate) 
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Frequency vCps; 

Fig.2.11 Displacement frequency response curve for stabilised 
soil (for circular plate) 















67 



500 540 580 600 

!Cps) 

abilfscd soil(for long strip) 


Circular Plate 



300 340 380 420 460 500 540 

Frequency (Cps) 

Fig.2.t5 Displacement frequency response curve for 
natural soil (for circular plate) 
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300 340 380 420 460 500 540 580 600 

Frequency (Cps) 

Fig.2.18 Displacement frequency response curve for natural 
soil (for long strip) 







m 


470 


For force 5.0 kg 

• • Natural soil 

X— "X Stabilized soil 
O — O Soil contained in bin 


1.0 

4.0 5.0 10 16 20 30 


L/B ratio 

Fig.2.19 

Variation of frequency with 


L/B ratio 





Maximum deflection (in mm) 


7.3 



L/B ratio 


Fig. 2.20 Variation of maximum deflec- 
tion with L/B ratio 



IL, 



.005 .009 .013 .017 .021 .025 .029 

Maximum displacement (in mm) 

Fig. 2.21 Variation of maximum displacement 
with force 



Natural soil 



Fig. 2.22 Variation of spring constant ratio with 
L/B ratio 
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O 



Fig. 2.23 Variation in Damping factor with L/B ratio 


6.85 kg/cm 



Fig* 2.24 Displacement frequency response curve 
for a beam resting on ideal springs 




Maximum displacement (in mm) Maximum displacement (in mm) 


Forcer 5.15 kg 
K =s 6.85 kg/cm 
I r 0-867 cm^ 
L = 50 cm 


. 

Theoretical 

Experimental 


Distance from origin^^ 

2.25 Maximum displacement variation 

of beam at different 
encies for force 5*15 kg. 
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50cps 
- ,1 
-.2 

- . 3 

4 

- -5 
-•6 

55cps 


- • 7 

BOcps 
- . 1 

- -2 
-•3 

- • 4 

- -5 

65 cps 


70cps 
- .1 
-•2 
- -S' ■' 
80 cps 


90 cps 

::V 


0.0 0.2 0.4 0.6 0.8 1.0 


Fig. 2.2 6 Ndxirnum displciccmeni va along the 

different 

frequencies for force 9.40 kg. 






Maximum displacement (in mm) Maximum displacement (in mm) 




GH^TEH 3 


DYIUMIC RESPOKSE OF BE&S ON FLiSTIC POUIfDiiTION CiREYING A 

CONCENTEiiSED M/SS 


5.1 GENEEiL 

This chapter deals with the analytical solution of Euler-Bernoulli 
heam, carrying a concentrated mass attached to it, resting on elastic 
fotmdation and subjected to a steady state dynanic force. This 
problem has a lot of bearing on practical situations. A machine 
resting on deep beam, a coltunn on strip footing, a mechanical 
oscillator (to produce a high force oscillation) on footing* for the 
determination of soil properties, are the few examples which simulate 
this problem. 

In the literature, solutions for beams carrying concentrated niass 
and simply supported at the ends without foundai;ion are available ' 
(12, 2, 15* 49» 22). So far no solution is reported for beams on 
elastic foundations with concentrated mass attached to it. In this 
chapter, the equations for displacement of beams with concentrated 
mass attached to it at the centre of the beam and resting on elastic 
foundation are developed using the Pasternaic foundation model and the . 
Winkler foundation model. The damping in the foundation is included 
for the analysis. The theory and the results are presented for beams 
of finite and infinite length. 

The general solution for any arbitrary position of concentrated 
mass on the beam can easily be developed by adding suitable boundary 
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conditions and continuity conditions in the derivation. However, 
solutions for "beatn with concentrated nass at the centre of the hesan, 
resting on Pasternak and Winkler foundation have been dealt here. 
The results are presented in non-dimensional form and the solutions 
are compared with existing solutions wherever possible. The 
theoretical results are also compared with experimental values 
obtained by conducting experiments on a very large channel , resting 
on elastic foundation of natural soil deposit and attached with a 
mechanical oscillator (of considerable mass) at its centre. The 
oscillator generates a frequency dependent harmonic force. 


3.2 GEmiL THEOEY 


The differential equation of transverse vibration of an Euler- 
Bemoulli beam on elastic foundation subjected to a dynamic load can 
be written as - 


El 


3^w 


- G. 


2 

3 w 


3X 


H. Kw + o If + pi 


at 


* q(x,t) 


(3.2.1) 


where 

E is Young's modulus of elasticity, 

I is moment of Inertia of the beam, 

G-^ is shear modulus of the soil medium, 

K is spring constant of the soil medium, 

c is daii5)ing coefficient of the soil medium, 

pA is mass density of the so il mediunt ; 
c 

w is verrbical displacement of beam, 

El is modulus of rigidity of beam, 
and q(x,t) is dynamic forcing function. 
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In equation (5.2.1 ) the soil medium is idealized as a Pasternak 
foundation model characterized hy the foundation constants and K. 


Defining the constant and the variahles as 


X » 




X 

X 


Sill 


0} 


47-^ ,z=7,y = “!-.r* 


X ’ 


m » e - t, 0 ^ = 2 pA^ 
n 


D = 


/M , 


(5.2.2) 


•e = mt , 


where 


a,nd 


X is the characteristic length, 

z is non-dimensional distance variable, 

y is non-dimensional Displacenent variable, 

r is dimensionless forced frequency, 

9 is non-dimensional time variable, 

c. is critical Dancing constant, 

D is damping factor, 

n is characteristic shear value of soil medium. 


Equation (5.2.1 ) can be written in the non-dimensional form as 


+ 4y + 8lbc ■^ + 4r^ ^ g(z»6) 


(5.2.5) 


For the unloaded surface of the foundation, outside the finite 
beam, q(z,e) will be zero. Thus the equation of motion of this portion 
can be written as 

2 

9 w « aw- 

So - K .£ + o 5^ - 0 (5.2.4) 

3x 

where w^ is the foundation deflection beyond the length of the beam. 
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Equation (5.2.4) can be e^cpressed in nondimensio ml form as 


3 

^ (1 + 2Dr) = 0 (3.2.5) 

Sz X 

where y^ is the nondimensional displacement of foundation beyond beam. 
The general solution of equation (5.2.5) can be obtained as 


yf(z» e) 



sin 2Dr ^ z 



sin 21b: ^ z 

... (5.2.6) 


ire 

e 


5.2,1 Boundary Conditions for Finite Seam 

Pig, (3.1b) shows a uniform finite beam, supported by a 
Pasternak foundation, a mass 'M' and a sinusoidal force P = P 

o 

is concentrated at the middle point. This point is taken as the 
origin of the coordinates. Because of symmetry, only the right side 
of the beam is considered. At the origin, the slope of the deflection 
curve is zero and the total shear force applied is equal to shear 
force restored. The shear force applied will consist of the steady 
state constant exciting force P^ e at this point, and the inertia 
force, due to the concentrated mass M attached at the point. The 
restoring shear force will be provided by the beam and the shear 
layer of the Pasternak foundation. In mathematical form, these can 
be written as 


at z = 0 


M. 

3z 


0 


and 


El 


.2 ■ 
32 ' ^ 






2 

X 3 : 0) 


2 

n 


y 


(3.2.7) 
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or |z “fll 4a R at z = 0 (3*2.8) 

3z 

At the free end of the beam, the bending moment will be zero, 
deflection of the beam will be equal to deflection of the foundation 
just away from the beam, and total shear force will be equal to the 
concentrated foundation pressure. In mathematical form these boundary 
conditions at the free end can be written as 


at z * 

= R 

1! 

O 



(5*2.9) 



3z 






y(R) = yf(R) 


(5*2.10) 

and 

El 

" x2 


MM 

M 

32 

o 9Z 



A 

_3z‘^ n 




or 

5 

T - 

32 

if = i? 

9z 

sy^ 

3S 

at z « E 

(5*2*11) 


At infinity, the deflection of the free surface of the 
foundation must vanish. Thus 

y(f)^ » = 0 (3*2.12) 

Using these conditions (5*2.7 to 5*2.12), the basic differential 
equation is solved by Laplace transform. 

5 .2 .2 General Solution for Finite Beam 

The general solution of Eq. (5*2.5) can be taken as 

i(re - ^) 


y(z, e) = y(z) e 


(3*2.13) 
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Putting this solution in basic differential equation for homogeneous 
solution, (5.2.3) yields 



2 2 

1 ^ 3 y „i(3::0-i})) , . i(re-’<t)) ^ „ -n- i(r6-A) 

4 ~2 — 2 ® +4ye ^^ + 8 Dir ye ^ 

n 3z 

+ 4('-r^) y - 0 (3*2.14 ) 


After arranging the equation (3.2.14), it can be written as 
422 

^ „ 4 J_Z + 4(1 _ / + 2 Dir)y = 0 (5.2.15) 

3z^ n 3z 

For maximum deflection in a time cycle 

« / (1 - /f + {2Drf ( 3 . 2 .I 6 ) 

Thus , equation (3. 2.15) can be expressed as 


4 2 2 

^ ^22 

3Z n 3Z 


+ 4y y = 0 


(5.2.17) 


For the solution, the Laplace transform of equation (3.2 .17) can be 
written as 
4 

L(^) = y(p) - P^ y(o) - p^ y*(o) - py^’(o) - y«"(o) 

3z 

2 p 

“ P Kp) - P y(0) - y'(0) (3*2.18) 

3z 

L(y) " y(p) 

where L(f) denotes the Laplace transform of the function F 
(p) is transform variable. 
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On putting these transform values in equation (5«2.17) and. solving 
for y(p), the resulting expression is 

(p^-A^p) y(o) (p^-^^)y'(o) py"(o) y"'(o) 

y(p) = — 7^ — ?; s — TT- + — s— 7? - + — s — =- + 


(p^+C^)(p^-3)^) 


a" a" 
(5.2.19) 


where A = 2 X /p 


„2 _ A^ + /7T4? 

" 2 


(5.2.20) 


I? = 

a Q- 

B = 

Inverting the Laplace transform, equation (5.2,19) can he written as 


y(z) =. “2” 5 cosh Lg^z + cos Cg^z)y(O) + q p (C^ sinh 


D z + 
a 


+ I? sin CgZ ) y'(0) + (cosh LgZ - cos Cgz) y"(0) 

+ cV (^a V - \ y"’(0)l 

a a -i 


(5.2.21) 


The boundary conditions can, now, be expressed as 
at z = 0 , = 0 

aZ 

p x^ 

y«H(0) - A^ y.(o) = .4a Hr^ y(0) - 


Thus equation (5.2,8) can be written as 

T, 2 

F X 

y»'»(0) « -4a Sr^ y(0) - ° 


(5.2.7) 

(5.2.8) 


2EI 


(5.2.22) 
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y”(R) » 0 

(3.2.9) 

y(R) “ yf(R) 

(5.2.10) 

y'’'(R) - ^ y’(R) = yf(R) 

(3.2.11) 


On solving these houndary conditions, the initial parameters can he 
obtained as 


yi'>(0) » -4a Er^ y(0) - 


(5.2.22) 


y’(o) 


^ sinh DR + C sin C R 

O d a» 8. Q 

I? cosh D„R + C? cos C R 


+y(o) 


i4a Rr 


(D sinh D R + C sin C R I?(cosh D R-cos C R 

^ ci 3i SL cL M & ct 

■ D^ cosh D R + cos C R cosh D R+C^ cos C R 

a aaaa aaa 


and 


y’(o) 




(3-2.25J 

(5-2.7) 


2EI 


y(o) = 


V “a® + 

airi D E ain C E}+ 2(I?+oh(D?-oh''/^ ( 

a, Ot Ei B, o> a, 

(D cosh D R sin O R - C sinh D R cos C R) 1 
^ a a a a a & ' JL 

oosh \S cos C 3 ^E(o|+I?) 2 1 1+{2Ilrf 


- 4<i Ex^ x{1+(2arf (I^ oosh B^S sin C^H 

- "a °aE) " 4« <2<=M * 

cosh Dg^R cos CgR + Cgl^(D^-C^) sinh D^R sin CgR} (5-2.24) 
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M the y( 0 ), y'( 0 )j y"( 0 ), y*"( 0 ) (the initial parameters) 
are known, the displacement at any point along the length of the 
beam can be determined using equation ( 5 . 2 . 21 ). 

5.2.3 Slope, Bending Moment, Shear Force and Soil Heaotion 
Expressions for Finite Beam 


The expressions for the slope of the deflection curve, the 
bending moment, shear force and soil reaction for the finite beam* 
on Pasternak foundation can be derived from the displacement equation 
(3.2.21). These are as follows. 


Slope of the deflection curve 


S(z) 




sinh D z - D sin C z) y( 0 ) + 

cL o 0,1 


C^a sinh Dg^z + sin G^z) y''( 0 ) + (cosh D^z-cos C^^z) 

y'"(o)J (3.2.25) 


Bending moment of beam 


'z) 


El 




C^D^(oosh D z-cos C z)y(O) + (l^ cosh B 

8tf Cl, d a S, 


+ sin Cgz)y"( 0 ) + (B^^ sinh sin Cg^z)y>"( 0 ) 


(5.2 .26) 


Shear force T(z) 


El 






+ (B^ sinh B z-G^ sin C z)y"( 0 ) + (iF cosh B, 

8. Q, 3- 3 t 


+ G cos C„z 
a a 


) y”(o)^ 


(5.2.27) 
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Soil Reaction pCz) = 


Si 3/ 


+ (D^ cosh D z-C^ cos C_z)y^'(0) 

0f Ci ^ 


+ (D^ sinh I)„z - sin C a 

S' 3* a> 


3)y'”(0) (3.2.28) 


3.5 FUTITE mM on TiTBTELER FOURRiE'ION 


The solutions for deflections, slope, bending nonent, shear 
force and soil reaction are presented for the finite beam on Winkler 
foundation. Making the shear coefficient term, A, zero in equation 
(3.2,1), the differential equation of transverse vibration of the 
beam, can be written as 


El + Ktt + 0 + p A^ = <l(x,t) 

3 x 3 *t? 


(5.3.1) 


And in non-dimensional form it can be expressed as 


+ 4y + sDr + iifT 
3Z^ 38 


^ + 4r^ ^ « |:= q(z,e) 

38 "1 


(5.5.2) 


For free-end beam, the bOLindary conditions can be listed as 


at z = 0, = 0 


(3.5.5) 


at z = 0 , = 

3 Z^ 


'■ " ■ -om ' T + 4 a Hr y 


(3.5.4) 


at z = R, ^ = 0 


(5.5.5) 


(3.3.6; 
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-Assuming the general solution as 
y(z5e) = y(z) 

and solving equation (5»3*2) for boundary conditions (3 *5 *3 "fco 5»5*6)» 
the expression for general solution can be obtained as 


1 1 

y(z) <= (cosh gz + cos gz) y(0) + -^ (sinh gz + sin gz) y*(0) 


1 i V 

+ — 2 (cosh gz - cos gz) y"(0) + — = (sinh gz - sin gz)y‘"(0) 
2 & 2g^ 

... ( 3 . 3 . 7 ) 


where 


F X 

yii.(O) = _4a Rr^ y(0) ~ -|gj- 


(3.3.8) 


y"(0) = y"'(0) 


's(r) u(r) - T"^(R) 

.s(r) t(r) - ij(r) v(r) J 


(3.3.9) 


yt(0) = 0 


(3.3.10) 


y(o) 


F^ X^ T(e) Y(r) - S^(e) 

[^g^{U(R) Y(r) - S(r) t(e) } + 4 a {S^(R) 


-t(r) 1{' R )}’2 


... (3.5.11) 


Thus the general expressions for displaceuBnt of the beam, slope of 
of the deflection curve, bending moEBnt, shear force and soil reaction 
can be expressed as 


[{Y(R) t(r)-s^(r)} s(gz) + {s(r)u(r)-t^(r)} n(gz) 

F^ X^ + g^ {S(r) T(r) - 13 ^ (R)> V(ez) ] 

{jP {S(r) T(r)-U(r)V(r)} + 4 a Rr^ { T(r)Y(r)-S^(r) } ]] 

(5.3.12) 
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Slope of the deflection curve 


S(z) = Y(ez)e y(0) + S(3z) yt(0) + i T(3z) y”(0) + -^ TJ(3z) y'"(0) 


(5.3.13) 


Bending moment 


El L2 


u(ez) y(0) + 3V(3z) y'(0) + SCgz) y”(0) + -1 T(3 z) y*'’(0) 

... (3.5.14) 


Sheer force 

T(z) “ 'l(ez) y(0) + 3^ U(3z) y'(0) +!3'^(pz) y"(0) + S(ez)y'"(0) 

... (3.3.15) 


Soil Reaction 


p(z) “ S(3z) y(0) + ^ 'r(ss) yf ( 0 ) + 3^ R(3z)y‘’(C)+ 

+ 3V(3z) y"'(0)1 


(5.3.16) 


■jfhers 


S(Sz) 


cosh gz + cos gz 


T(3 z) 


ij(ez) 


V(3 z) = 


sinh 3z + sin 3z 


cosh 3z - cos 3z 


sihh 3z - sin 3Z 


(3.3.17) 
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3.4 INFINITE BEAM ON PASTSBUiK FOMDATIOI? 

The solution for infinite beam can be obtained from the differential 
equation (5.2.5) by changing the boundary conditions suitably. The 
solution can be expressed as 


L z 

-D z 

D z 

sl 

y(z) = e cos C^^z + 0^ 

a 

e 

e-D Z 

a 

Sl 

cos Cg^z + e sin C^z 

+ c. 

e 

sin C„z 

4 


a 


C„, are defined by equation (5.2.20). 

3r' cL 


For infinite beam at z ■*■<*> , the displacement, slope, Bending 

moment, shear force will have some finite values, thus the terms 

D 

z 

containing e is to be ommitted. At z = 0, the boundary conditions 
will be as follows 


and 


at z = 0 , = 0 

92 


9z' 


9z 


4a r y - 


2 EI 


( 54 . 2 ) 


( 3 . 4 . 5 ) 


Using these boundary conditions in equation ( 3 . 4 .Oj the solution 
can be obtained as 


y(z) = - 


^o ^ 


1 r “V ^ 

o-PT 9 — 5 5 ® sin C z 


& 


G 


... (5 .4 .4 ) 


and expressions for other functions can be written as 

-D z 


slope S(z) 




sin C z 

3. 


2EI 




( 54 . 5 ) 
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Bending Moment M(z) = - 




(C^ + I^) 
^ a 


^ (SVX + 2a Aj 


^ "®a ® ^a® ^a* } (5*4*6) 

(o' + 1?) 


-r z 


Shear force T(z) 


a a' 


}0a\(0a 4) + 2a ^2 cf 


-2 2 “V «■ 

(ih - C )e sin C z - 2C„I)„ e cos C„z 

da, '' 


-D.z 

a'a ^ — a’ 


and Soil Reaction P(z) 


2X 


X (C^ + !^) 


20a“a(ot+4) + 2« 0 


ic 


a 


r-(^ 

^ a 


J 

... (5.4.7) 


-D z 

„ e sin C z 
a a 


■z 

+ e cos C^z) + 5CaPg^(Cg_ sin C^^z+Bg^ cos C^z) 


... (54.8) 


5.5 BIRIHITE BEAM ON WINKLER FOIMDiflllOH 


If the shear coefficient is made zero, then equation (5.4.4) 
can he modified to give the deflection equation for Winkler foundation. 
Thus 


y(z) = 

0 

-8z / . \ 

e (,cos gz + sin gz; 

2EI 

(48^ + ar^) 

S(z ) = 

^0 

^ « ~6z 

2 g e sin gz 

2EI 

(4g^ + ar^) 

M(z) « 

P X 

2g e ^ (cos gz - sin gz) 

~ 2 

(4B^ + ar^) 


(5.5.1) 


(5.5.2) 

(5.5.5) 
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(e*®^ oos gzj 

^(z) = +— , 3 _ ^ 

(4& + ar ) 


P(z) = - 


P e Tcos Qz + sin gzl 


( 3 . 5 * 4 ) 


(3.5.5) 


5.6 RESULTS 

The response curves are plotted for displacement, frequency etc » 
and are presented in dimensionless form. The equations derived for 
dis-nlacement, slope, Lending moment, shear force and soil reaction fcm 
these two foundation models can be used for the response curves. Various 
results of the responses are presented below. 


3 .6.1 Finite Beam on Winkler Foundation 

The response curves are plotted in terms of displaoenent ratio 
and frequency ratio. Displacement ratio can be defined as the ratio of 
displacenent at any frequency ratio, divided by the displacement at zero 
frequency ratio. Figs, 3,2 to 3,7 present the displacement-frequency 
response of finite beam of various lengths, for various values of dancing 
factor and concentrated mass. Some of these figures present the 
resnonses at various points along the length of the beam. 

Fig. 5.2 presents a typical response of a finite beam of dimensionless 
length S « 0 , 5 » with no concentrated mass attached to it. The response 
cuirve is for z » 0,3 (dimensionless distance from origin x/jl). The 
variation of disnlacement with frequency is plotted for various values 
of damping factor D (O.O, 0,05, 0,10, O.15). While the maximum displacement 
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ratio is infinite at resonant frequency for undamped foundation, the 
mA-v -iTrmm displacement becomes finite and decreases with increase in 
dan^jing factor. The response at higher frequency, approaches the 
same magnitude of displacement for any damping factor. 

The response of a beam of E = 0.8 with concentrated mass is 
quite similar to that in the Eig, 5 . 2 . The response curves shown in 
Fig. 3.3 are for z = 0 and R « 0.8. The displacement ratio is infinite 
for undamped foundation and no concentrated mass attached to the beam. 
However, when the concentrated mass is attached to the beam, the maximum 
displacement becomes finite even without damping. The maximum displacement 
is decreasing with increasing ot (concentrated mass ratio). M in the 
Fig, 3 , 2 , the displacement values tend to approach the same value at 
higher frequency. 

For beam of unit nondimensional length (R = 1.0), the displacement 
at the centre (the point at which the force is applied) of the beam with 
frequency ratio is plotted in Fig. 34 . For various values of a (from 0.0 
to 0 . 6 ), the displacement-frequenjoy response are plotted, talcing the 
value of damping factor equal to 0,1. Maximum displacement can be 
observed to be decreasing with concentrated mass ratio, just as in Fig. 3 ♦3* 
The displacement values approach the same value at frequency ratio 2 . 
if ter that, the displacement decreases fast with hi^er concentrated mass 

ratios. 

Figs. 5 , 5 , 3 . 6 , 3 . 7 , present the displacement-frequency response 
for beams of nondimensional lengths R *» 0.6, 0.8, 1.0 respectively. These 
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curves are for the case in which the foundation has the damping factor 
of 0,1 and the value of concentrated mass ratio a 0,3. Each figure 
shows the displacement-frequency ciorves at many points along the length 
of the heam. Maximum displacement occurs at z = 0 and with the increase 
in z, the maximum displacement decreases, At one particular frequency 
ratio of 1,4 » "the displacement along the length of the "beam can 
he seen to be the same, whatever he the length of the heam. Beyond 
this frequency, the displacement increases with z. While the displacement 
at the centre of the heam decreases with increase in frequency, the 
displacement at other points on the beam, increases slowly with the 
increase in frequency. With increase in z, the rate of increase of 
disnlacement increases with frequency. For beams of R = 0,8 and 1.0, 
the displacement values, for frequency ratios between 0 to 1.0 are less 
than the displacement values at frequency ratio 0, For z less than 0,8, 
the disnlaoement values in this frequency range are always greater than 
the value at frequency ratio 0.0, For R = 0,8 and above, the value 
of disnlaoement can be seen to increase beyond the resonant frequency 
ratio. 


From these basic response curves, the maximum displacement 
variation with damping can be plotted for various concentrated mass 
ratios. Fig. 3*8 shows the variation of maximum displacement with 
damping, of a finite beam (R = 0,6) at its centre for various values 
of concentrated mass ratio (a= 0.0, 0,1, 0.2, 0,3). It can be 
observed that, with increasing damping, the maximum displacement is 
decreasing. The rate of decrease is high for the beam with no 
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concentrated mass attached to it. This rate comes down with a increasing. 

Fig, 5,9 shows the maximum displacement variation along the length 
of the beam resting on Winkler foundation v/ith no dancing for various 
values of concentrated mass ratios (a = 0.2, 0.3,0. 5, 0.6, O.7). -As 
the 2 increases, the maximum displacement decreases. This decrease 
is not proportional. For lower values of z, the rate of decrease is 
small. For the higher values of z, the rate is high. With higher values 
of a , the rate decrease of maximum displacement is small. 

This can be clearly observed from Fig. 3*H* This figure presents 
the variation of maximum displacement with concentrated mass ratio a . 

The beam length is E = 0,8, and foundation has no dan^jing. The maximum 
displacement variation at points along the length of the beam has, more 
or less, the same rate of decrement. 

For beams of finite length (E = 0,5, 0,6, 0,8, 1.0), resting on 
Winkler foundation and damping included, (D = 0.1 ), the variation of 
maximum displacement along the length of the beam for various values of 
concentrated mass ratio., are shown in Figs, 5»11, 5*12, 5*15* 

It can be observed from these figures that the maximum displacement is 
decreasing with z increasing. However, the rate of decrease is less in 
this damped case compared to the undamped case. With increasing 
concentrated mass ratio a , the rate of decrease of maximum displacement 
can be seen to be the same. 

Figs, 3«15» 5*17» 5-18 present the maximum displaceuent 

variation with concentrated mass ratio for finite beam of lengths 
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R = 0,5, 0.6, 0.8, 1.0 respectively, resting on v/inkler foundation with 
damping (D = .l). The decrease in maximum displacement with concentrated 
mass ratio is not proportional. M the value of a increases, the rate 
of decrease becomes small. For different values of z, the rate of 
decrease of maximum displacement is nearly the same. 

It can, again, be observed from these figures that the maximum 
disnlacement is increasing as the length of the beam (r) is increasing. 
Pig. 5*19 presents the variation of maximum displacement at the centre 
with the dimensionless length of the beam. These curves are plotted for 
various values of a(0,1 ,0,2,0,3)* Prom the figure, it can be observed 
that the increase in the maximum displacement is very marginal. 

3.6.2 Infinite Beam on Winkler Poundation 

Por infinite beam on Winkler foundation with damping (D = 0*1 ), 
a typical displacement frequency response for a ® 0,1 is shown in 
Pig. 3 *20, for various values of z. The maximum displacement (not the 
ratio) occuis at the centre of the beam at frequency ratio 1, This 
maximum displacement decreases along the length of the beam. With 
increase in concentrated mass ratio, this maximum displacement, again, 
decreases. Pig. 3*21 shows the variation of maximum displacement 
with concentrated mass ratio for damped and xmdamped foundation model. 

In undamped foundation model, the displacement at resonant frequency 
is quite high, thus the displacement corresponding to frequency ratio 
0.9 is presented in this figure. The decrease in maximum displacement 
is almost linear. Pig, 3*22 presents the maximum displacement variation 
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for damped and undamped foundation model of the beam with z values* 

Out of the two curves, one curve is for the beam not carrying the 
concentrated mass and the other for beam with concentrated mass attached 
to it* (a SI 0*1 )* is usual, the maximum displacement decreases slowly 
at lower values of z and then faster at higher values of • 

5*6*5 Infinite Beam on Pasternak Foundation 

Using the equation (5*5*12) for infinite beam the response curves 
for deflection of beam resting on Pasternak foundation are sho?m in 
Pigs. 5*25 to 5*28 for various cases* Displacement-frequency response 
curves are presented for varying values of concentrated mass ratio, 
shear coefficient A of shear layer of Pasternak foundation, damped and 
undamped cases and at various points along the length of the beam a^way 
from the centre# 

Pig. 5*25 presents displacement-frequency response of an infinite 
beam with no concentrated mass on it for damped and undamped cases, at 
points on the beam away from the centre# The value of shear coefficient 
A is 1 *0* Maximum displacements occur at resonant frequency and is 
maximum at the centre of the beam# As z increases, the displacement 
decreases# While in undamped case, the maximum displacement is 
reasonably high, in damped case, the m^imum value is very lev/# 

The displacement ratio is plotted against frequency ratio to present 
the response curve of an infinite beam resting on Pasternak foundation 
(As= 0#5, 0#l,aa 0 . 0 ) in Pig# 5*24* These response curves are 
for various values of z from the origin# A^ can be observed from the 
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figure, the maximum displacemeniJ-' ratio is quite high. The maximum 
displacement values decrease with increasing a. But this decrement 
is very small. At other frequency ratios, the displacement values "are 
very close to each other for all the v'lues of z. Only at resonance, 
there is difference in displacement values. 

Figs. 5»25> 3*26, 5*27 present the response curves between displacement 
ratio and frequency ratio for a = 0.00, 0,2, 0,4 respectively, with B = 0.0. 
The response curves have been shown for various values of shear coefficient 
A. For lower values of shear coefficient, the magnitude of maximum 
displacement can be seen to be very high. At higher values of shear 
coefficient, the maximum displacement ratio is low. It can also be 
observed from these figures, that at frequency ratio 14 , the displacement 
ratio is same for any value of shear coefficient A, At higher frequency 
ratios, the displacement ratio values increase with increasing shear 
coefficient. 

Fig. 5,2R shows a typical set of response curves for foundation- 
having damning ratio 0,1. With shear coefficient increasing, maximum 
displacement ra-tio decreases very fast. It can again be observed that 
the displacement ratio at the frequency ratio 1,4 is the same for any 
value of shear coefficient as in the undamped case. The displacement 
ratio values are higher for higher values of shear coefficient at higher 
frequencies. 

Variation of maximum displacenent ratios with concentrated m^s 
ratios for various values of shear coefficients axe plotted in Fig, 5*29 
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for D = 0.1, It can be observed, that the maximum displacement ratio 
decreases ^ith higher values of concentrated ;.sss. But as the shear 
coefficient becomes high, the effect of hi^er concentrated mass on 
maximum displacement is less. For 1.0, the maximum displacement 
decreases reasonably with concentrated mass ratio. For A - 5»0, there 
is hardly any decrease in maximum displacement ratio with a . 

Maximum displacement along the length of the beam decreases, as 
can be seen from Fig, Response curves are plotted for 

various values of shear coefficients, The beam is attached with a 
concentrated mass of value a = 0,1, Foundation has the damping D = 0.1, 
The rate of decrease of maximum displacement for z values is nearly sane> 
but the reduction in the maximum displacenEnt is quite high with higher 
shear coefficient.' 

In figure 3*31> curves are drawn showing variation of maximum 
displacement of the centre point of the beam with shear coefficient A. 
These curves include the cases of undanped and damped foundation and 
a «= 0,0, and 0,2, The maximum displacement decreases very fast for 
undamped and damped cases. The maximum displacement for undamped case 
is higher than for damped case, ibove shear coefficient value A= 2.0, 
the maximum displasement for damped and undamped case become almost equal. 
The effect of concentrated mass is appreciable at lower values of shear 
coefficient, 

3 ,6 ,4 Finite Beam on Pasternak Foundation 

The response curves of a finite beam on Pasternak foundation 
can be drawn, using the deflection equations (5. 2 ,24 )• Figs, 3*52 and 3*53 
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present the displacement frequency response curves for various values of 
concentrated mass ratio. The heam is free at the ends and the concentrated 
mass is attached at centre of the beam. A sinusoidal force is applied 
at this point. The nondimensional length of the beam is O.707I • The 
results are plotted for damped and undamped cases. 

Fig, 3.32, presents the displacement-frequency response ciirves for 
damned and undamped cases for beam resting on Pasternak foundation. 

For A ts 0, the foundation will reduce to Winkler model. The response 
for A = 0 is compared with the response for A = 1 ,4'l41*^hile for A = 0, 
the maximum deflection occurs at frequency ratio 1.0> for A= 1,4141, the 
resonant frequency shifts from 1 to 1,3 • Rades' solution gives the same 
response for R = 0,7071 and A *= 1.4141, 

In Fig, 5*35> Displacement response cuives are drawn for shear 
coefficient A 1.4141, damping factor D = 0,1. These displacement 
values are at the centre of the beam, Ag in Rig. 3»52» the maximum 
displacenBnt occurs at frequency ratio of 1,3. With increasing concentrated 
mass ratio, the maximum displacement decreases. At higher frequency 
ratios, the values of displacement decrease and for any value of a , 
these are nearly equal. 

The responses were calculated for shear coefficient values of 0,5 
and 0,7<171 also. The variation of maximum displacement is shov/n in 
Pig. 5»54 f or ot ■ 0. 0,0,1, 0,2 for various values of shear coefficient A. 

The rate of decrease in maximum displacement of the beam with no 
concentrated mass is small at lower shear coefficient values. At higher 
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values of 1, the rate increases. the concentrated mass on the beam 
increases, the rate of decrease tends to he uniform. 

Pig, 3 •55 shows the variation of maxinaim displacement along the 
length of the beam for P = 0 , 7071 , a = 0,1 and D = 0 , 1 . vVith ’z' 
increasing, the maximum amplitude can be seen to be decreasing. The 
rate of decrease for different values of A seems to be similar. 

The general trend in the variation of maximum displacement with 
other variables for finite beam is similar as was observed for the 
infinite beam on Pasternak foundation. 

5.7 YERIFIC ATIOIT OF SOIIJTIONS MI) BESULTS 

In this section, the general solution for beam carrying concentrated 
mass at its centre and resting on elastic foundation is compared T<ith 
the existing solutions for beam, resting on Pasternak foundation and 
beams carrying concentrated mass, simply supported at the ends. The 
analytical results of the solutions are compared with experimental 
values, obtained by conducting the experiments on a long beam with 
concentrated mass attached to it and resting on the soil medium, 

3 *7 *1 Verification of General Solutions 

The displacement expression for beams on Pasternak foundation 
and Mnkler foundation are given by equations (5,2.24), (5.5.II). If the 
concentrated mass term a is made zero, these expression will be as follows. 
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The displaceraent equation for team on Pasternak foundation is, 

p r 

- -Iir V <=“= 

sii* D^K ,iaC^R + 2 

(l + (2Dr)^ cosh D^E sin C^R-C^ sinh B^R cos 
c\ D^(C^+B^)(b! cosh B R sin C R + sinh B R cos G E) ■ 

3. 3 rt, 3 3i 3. 3. 3. 3. 3. 


'”c 


cosh B^ cos C^R(cJ+4f 2(B2-c2)''/2(^+(2Brf 


... ( 3 . 7 . 1 ) 


and the displacenBnt equation for "beam on Vifinkler foundation is 


(sinh gR + sin BR) (sinh gR - sin BE) - ^ (cosh gR 


y(o) = 


Pq X^ + cos gE)^__ 


|”g^ (cosh BR - cos BE) (sinh gR - sin BR) - 
^ (cosh BE + cos gR) (sinh gR + sin BE) } "^ 


... ( 5 . 7 . 2 ) 


The same (iq. (3*7.2)) can he simplified as 

-P 


y(0) o - ° Ql+cosh gR cos gR^]/ Q’osh gR sin gR + sinh gR cos BR^ 

2BI g^ 

... ( 3 . 7 . 5 ) 


Rades ( 40 ) has presented the solution of Euler-Bernoulli heam of finite 
length on Pasternak foundation and subjected to steady state dynamic 
force. The deflection expression for beam is given by 
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-4 , r4 


+ X"^) cosh X cos V + p X (x^ - Ij"^) Sinh x sin y 


-2' 


y(o) 


y D -i.jiJ.iA A 


where 


2 —2 —2 —4 —4 — i - — 

■F^L + 2v X } - Y(y - X )(X cosh X sin y 

• 2gT^t^2 —2 —2 —2 —3 — — -*3 — — > 

, y X (y + X )(X cosh \ sin y + y"^ sinhx cos y, 

«. «,2 —2 — 2—2 2 — — 

- y (X - y )(y + X ) cosh X cos y 


^ .w 


(3.7.4) 


F^ = Force applied 
l - (2iif + 2(m‘^ + 

- (.2i + 2(n^ + n‘*)'/2)l/2 


m 


('£if.xi/2 

'‘4EI ^ 


Y = 


n 




For Winkler foundation, the expression is given by 
2 

-F^ L (cosh 5 cos 5 + 1 ) 

° 2EI ^ (cosh 5 sin 5 + sinh 5 cos 5) 


(5.7.5) 


Term in the expression (3.7,2) are corresponding to y , x respectively. 

Damping is zero in Rades equation and y ocrresponds to 

T - — f— T-r/ p ( 3 - 7 - 6 ) 

Equation of beam on Winkler foundation (3 *7 .5) is same as equation 
(5.7.5). ^he term 5 is equivalent to gE. Thus the general solution 
of a finite beam carrying a concentrated mass on Pasternak foundation 
boils down to Rades equation (40) if the concentrated mass ratio a is made 


zero 
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The general solution of the beam with concentrated mass can 
be derived for other boundary conditions. The solution of the beam 
deflection and frequency of simply supported beam carrying a concentrated 
mass at its centre is presented by Das and Deshmukh (15)» 'the derivation, 
the effect of rotary inertia and shear of the beam are also taken into 
account. If the rotary enertia and shear effect is neglected, the 
freauency equation can be written as 


1 Tb .... De - 
m 2 


(g tan - 6 tanh •^) 


(3.7.7) 


where M is concentrated mass attached to the beam and 


b = 




c ^2 2 
L 10 




Simplification of Eq. (3*7.7) gives 

(tan ^ - tanh ~) j " - jf =4 (3.7.8) 

The general solution for beam on Winkler foundation is given by 
equation ( 3 . 3 . 7 ). For boundary conditions of sinply supported case, 
this solution can be used to derive frequency equation. 


Boundary conditions are 


at z =» 0 ; *1^ = 0 , = 4a Hr y 


and atz = R*y»0» 

0 


2 

s y 


3z 


(3.7.9) 



The frequency equation can he written as 


y(0) (cosh 6R + cos eR)(sinh 3 R + sin gR) + y"(0) (cosh gR - cos gR) 

26 

(sinh bR + sin 3 R) + (sinh^ gR - sin^ 6R) = 0 (5»7-lO) 

26 

If there is no foundation helow the beam, then 6 = ^ . Solving equation 
(3»7*10) for boundary conditions (3 •7*9) yields, 

2 ^ 

Mw (tanh rR - tan rR) + 4^1 r'^ = 0 

2 

or (tan rR - tanh rR) — ? = 4 (3*7»'l'' ) 

4EI r^ 

where rR = ^ , and — ® ' a ^r ' ^ 

^ 4EIr^ c 

Equation (3»7*1 'I) and (3»7*’^) can be seen to be identical. 

3 •7*2 Verification of Results with Experimental Values 

Field experiments have been conducted to verify the analytical 
results. -A- line diagram showing the experimental set up is shown in 
Pig. 5 . 56 . 

3 . 7 *2.1 Experimental set up 

Fig. 3 *56 (a) shows the line diagram of a channel section, lO.O m 
long resting on the natural soil deposit. The mechanical oscillator 
(Saraswati SEA 35^*) was fixed at the centre of the beam. This mechanical 
oscillator, which acted as a concentrated mass attached at the centre, was 
producing frequency dependent sinusoidal force. A variable speed motor 
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was connected through a flexible shaft to drive the oscillator. Fig. 3*56(b) 
the block diagram of test set up. The variable speed motor runs the 
oscillator. The nsgnitude of the force level was controlled by changing 
the eccentricity of the rotating mass on the shaft of the oscillator and 
the frequency was varied by changing the speed of variable speed motor. 

The an^litudes of maximum djrnamic displacement at the centre and the 
points along the length of the beam were measured through a system of 
vibration Pickups' (iffl 126,127) 3ud vibrationmeter (iffi 75?) • 
oscilloscope was used for the visual check of the forcing function. The 
Pickups used in this experiments were the velocity type. In this, the 
voltage induced in the coil is directly proportional to the i^elative 
velocity between the coil and the magnetic field. The signal from this 
pickup was fed into vibrationmeter which is essentially a voltage 
meter calibrated to indicate the voltage in terms of the velocity of 
vibration system. The displacement and acceleration could be measured 
as the vibration meter had an integrating circuit and a differentiating 
circuit. A brief description of the oscilloscope is given in Chapter 2. 

In this test, the mechanical oscillator was tightly fixed to the beam 
(channel section). This oscillator utilizes the centrifugal force of 
unbalaiKsed eccentric masses to generate a harmDnic force. Two shafts 
with equal eccentric mass and rotating in opposite direction and tirhen 
in phase, prodtice a vertical sinusoidal force. With the increase in 
speed of rotation the magnitude of force also increases. 

Hie oscillator which has the total weight 75 ‘4 was attached 
with beam at the centre, firmly. The beam data are shown in Pig. 5*36. 
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2 

The spring constant of the soil was 525 Isg/cm and the concentrated mass 
ratio a= 0,54. 

The tests v;ere run for tharee different phase angles of eccentric 
mass ( <(> = 40®, 55“» 70®). The speed was varied from 600 rpm to I4OO rpm. 
The displacement and the velocity were measured at the centre and points 
along the length of "beam at 1 m interval. For each eccentricity, the 
tests were repeated a number of times and average of readings were noted. 
The values of velocity were recorded for cross check. 

5 •? *2 .2 Results and Discussion 

Pig. 3*37 presents the magnitude of displacements with 
increase in frequency for various phase angles (4^®j 55* > 70® )» a-'t the 
centre of the heam. While the rate of increase of magnitude of displacemj 
is slow for phase angles 40®>55'’ the increase for "JO* phase angle is quite 
rapid. The variation of maximum displacement along the length of the heam 
is presented in Fig, 5*58, % in the previous figure, curves are drawn 

for 3 phase angles of eccentric masses ((j) = 40®> 55®* 70°)* Beyond 
z => 0,R, the displacements have been observed to be negligible. It can 
be observed from the figure that maximum displacement decreases very fast 
with distance. For 4)= 70®, the amplitude of deflection is very high 
compared to other values of phase angles ((j) = 40“ » 55“)* 

Fig, 5*39» presents the theoretical results for the beam 
resting on the soil medium, as described above. The curves are drawn 
for a = 0 and a = 0.54» and force levels of 274 % s-nd 210 kg. While the 
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experimental curves are for frequency dependent force, the theoretical 
responses are for constant force excitation. Solid lines present the 
theoretical curve for a = 0»54> the dotted lines are experimental 
curves. For experimental curves, the force is frequency dependent. 

Thus at points (A) and (b) in Fig, 3»39> the force level of theoretical 
and experimental results are same. From these points, it can he 
observed that the experimental values of displacenents are quite close 
to theoretical values. 
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Fig. 3.1 (a) A beam element on Pasternak 
'foundation 
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Fig. 3.1(b) A uniform beam on Pasternak 
foundation 
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Fig. 3.2 Displacement frequency response curves 
for beam on Winkler foundation 
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Fig. 3.3 Displacement response curve for beam with 

concentrated mass at th® centre of beam resting on 
Winkler founds 
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Fig. 3.4 Displacement frequency response for beam 
with concentrated mass and resting on 
Winkler foundation (R= 10)^^ 


Beam of finite length 
on Winkler foundation 
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Fig. 3.5 Displacement frequency responses at various 
points along the length of beam (RsO-S) 
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Fig. 3.6 Displacement frequency response at various 
points along the length of beam (R=: 0.8) 
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Fig. 3.7 Displacement frequency response at 
various points along the length of 
beam (Rs 1.0) 
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Fig.3*8 Variation of maximum dispfa- Fig. 3-9 Displacement variation along the length 
cement with damping factor of the beam for undamped case, 

for various cocentrated 
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Fig.3.10 Displacement variation along the Fig. 3.11 Displacement variation along 

length of beam for damped the length of beam for 

case (R = 0.5). damped case (R= 0.6). 
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Fig.3.12 Displacement variation along the Fig. 3.13 Displacement variation along the length of 

length of beam for damped case(R=0.8) beam for damped case (R = 1.0) 
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Fig. 3.14 Variation of displacement with 
concentrated mass ratio tor 
::undamped-'xase-::^-- , 
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Fig. 3.17 Var iation of displacement with 
concentrated rnass ratio for dam- 
ped case (R=0.8) V 
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Fig. 3 d9 Displacement variation with Fig*3.l8 Variation of displacement 

non-dimensional length R. with concentrated mass 

ratio for damped case (R=1.0) 
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09*3.20 Displacement frequency response of an infinite beam 
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Fig 3.21 Variation of maximum displacement Fig. 3.22 Variation of maximum dis- 

with concentrated mass ratio. placement along the length 

of th <2 boam . 
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Fig » 3.23 Displacement frequency response curve for an infinite 
beam on Pasternak foundation (variation along the 
length of the beam) A = 1-0 
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Fig. 3.24 Displacemeht frequency response for an infinite 
beai^ Pasternak foundation (A=0.5).(yaria1:lon 
along the length of the beam) 
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Fig. 3.25 Displacement frequency response for an 
infinite beam with no concentrated mass 
(variation with shear coefficient A) 
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Fig ,3 .27 Dlsplacennent frequency response of an 
infinite beam with concentrated m 
^ ^ (variation with the shear 

coefficient A) 
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Displacement frequency response for an 
infinite beam on Pasternak foundation 
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Fig. 3*31 Maximum displacement variation 
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Fig. 3.32 Displacement frequency response of a 
finite beam on Pasternak foundation 
(variation of concentrated mass) 
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6.0 


5.O1 


Finite beam on Pasternak 
foundation 

R = 0.7071 
D= 0.1 

z = o.c 

As 1.4141 



Fig. 3.33 Displacement frequency response of a finite 
beam on Pasternak foundation (As1.4U1) 
(variation of concentrated mass). 
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Channel particulars 
weights 22.1 kg/m 

hs 20.0 cm 
b s 7. 5 c m 
Is 1819.3 cm^ 

As28.21 cm^ 

/long beam(channcl 
\ section) 



^Channel 

section 


mgW sin wt 


\ ^Mechanical oscillator 



Frame 

pickup 


^To vibration 
meter 


Natural soil 




Concentrated mass attached at the centre 
of beam 


Fig. 3.36(a) A line diagram of test arrangement 
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variable speed meter 
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shaft 
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Mechanical oscill 
ator 

Beam 

_J ^ 



jM. 


vibration 
meter. 
VELOCITY TYPE 


Cathode 

ray 

oscillo- 

scope 


VIBRATION 
TRANSDUCER 


TiTsTTfs 


Soil media 


Fig. 3. 36 (b) A block diagram for the test of a 
beam carrying a concentrated 
mass- at its centre . 
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CEAITER 4 


lENAMIC RESPOHSE OE BEAI® OR ELASTIC EOUmilOR 
SUBJECTED TO IMHJDSE LOADS 


4.1 &EEBRAD. 


The beams on elastic foundations subjected to impolse loads hare 
been the subject of great interest. The dynamic loading effects on 
stationaiy structures such as railway tracks, landing of aircrafts on 
runways, transient loads on pavements, footings subjected to impulse 
loads are some examples. In this chapter, the beam resting on. elastic 
founda'cions subjected to impulse loads is analysed, treating the soil 
medium as Ihsternak foundation model and Winkler foundation model with 
damping. 


The governing differential equation for such a problem taking icdio 
the consideration the so called internal damping and ejcternal damping 
has been given by Hoppoan (2l), as 


EI^ + 
8x 




5 

3 w , 3w 
3x 3t 


+ Kw + p A 


2 

3 w 
3t^ 


E(x,t) 


(4-1.1) 


where c 

e 

El 

I 

K 


is external damping coefficient, 
is internal damping coefficient, 
is modulus of rigidity of beam, 
is moment of Inertia, 
is spring constant of foundation. 




is unit mass of beam. 
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The second texm in the equation ( 4 . 1.l) is due to ejsternal damping, 
Hoppaan ( 21 ) again developed an expression for impi].se load in terns of 
coefficient of restitution e^, Por the impact of a sphere on beam, 
the contact force is given by 


f{t) = 



. irt 
sin — 

X 


for 


0 < t < 


(4.1.2) 


vitiere 


m is the mass of the sphere, 
s 

Vp is velocity of sphere just before striking the beam, 

e is coefficient of restitution 
o 

T_ is time of contact 

L 


Stadler and Shreeves (50) developed the differential equation with 
some modifications. Instead of external damping, the foundation 
consisted of a shear layer, besides the spring constant. !Phe differ- 
ential ecjiation is similar to the one described by Ihsternak (29) in 
his foundation model of soil medium and is given by 


El 


4 

3 w 
3x^ 


+ & 


^2 
3 w 

2 

3X 


+ Kw + c ■!“ + pA^ = P(x,t) 


2 

3 w 


(5.2.1) 


3t 


Stadler and Shreeves (50) gave the solution for impulse load at origin, 
as the general solution is difficult to be obtained due to complexity 
in deiiving the solution. However, by numerical analysis, the above 
differential equation can be solved for responses at any aibitraiy pjoint 
of the beam. 

In this chapter, the solution of the governing differential equation 
is attempted by Laplace transform* Several Laplace transform inversion 
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techniques were tiled for the solution, and Berger’s inrerstion 
tedinique (6) was used successfully for the general solution. Ihe 
results are presented for infinite and finite beams resting on 
elastic foundation. An approximate solution for an infinite beam 
resting on Winkler foundation with damping, using Post & Widder’s 
formula ( 4 ), is deiived and i«sults are presented. Solutions have been 
compared with existing results, 

4.2 G-EIIERAl BOBMJlAriOIirS . 

The basic differential equation of a Euler^Bernoulli beam resting 
on Ihsternak foundation is given by equation (3,2,3) and is given by 

- t? + 4y + SBr + 4 r^ -S-Z- = o (3.2.3) 

3 / 3 / 362 

!IIhe impulse load is applied at the origin at time t=0. This 
impiLse can be incorporated in boundary conditions with appropriate 
initial conditions at the origin. . He nee, applying the Laplace transform 
with jrespect to 9 , the equation (3.2, 3 ) can be written as 

4— 2— 

^ - A^^ + 4(1 + 2Dp + p^) y = 0 ( 4 . 2 , 1 ) 

3z 3Z 

wiere 

y( 2 ,p) = J y(z,0) e^®de ( 4 , 2 , 2 ) 

o 

p is Laplace transform variable 

y is Laplace transform of the displacement 
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4.2.1 Initial and Boundaiy Conditions 

(a) Initial Conditions . J\ist before the application of impulse 
load, the beam has zero displacement and slope of the deflection curve, 
lEhese initial conditions can be written as 


y(z,e) =0, at 6=0, ^ = 0 at 6=0 (4.2.3) 

O 0 

(b) Boundaiy Conditions for Unite Beam on Pasternak Foundation. 


The origin is taken at the centre of the beam. As the point 
is at the centre of the beam, the slope of the deflection curve is zero. 
The ^ear forces acting at this point are from the impulse load, ^ear 
resistance of the beam and the foundation. At the free eMs, the 
moment is zero while the shear resistance of the beam is balanced by 
the shear resistance of the foundation. In mathematical form, these 
boundary conditions can be expressed as 


at » = 0 , H =. 0 , am 

3z5 


2EI 


at z = R , ^ 

3z 3z 


Lx. ^ 


2»yf 


(4.2.4) 


9 z “ 9 z 

For unloaded surface of the foundation, the boundaiy conditions are 


3z^ 


(l+2Iir)y = 0 and as z 


= 0 


(4.2.5) 


(c) Infinite Beam on Ihsternak Foundation . The boundaiy conditions 
at the origin are the same as for finite beam. At infinity, the response 
functions have finite values. These boundaiy conditions, in mathematical 


form, can be written as 
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at z = 0, -^=0, and ^ = Imp, 

» 82 ’ 3^3 8 2 


9 

X 0) 


n 


2EI 


2 3 

and as z -*■ “ , ^'*^”’"8^ ^ ^~2’ ^ ^ 

82 3z 


(4.2.6) 


(d) jOliiLte Beam on Winkler Eoundation. For Winkler model the 
shear coefficient A is zero. Accordingly, the boundary conditions 
can be written as 


at z =: 0, = 0, and ^ = Imp 

8 z 

and 

2 3 

at 2 = R , = 0 , and = 0 

8z 8z 


( 4 . 2 . 7 ) 


(e) Infinite Beam on Winkler Foundation. The boundary conditions 
at z ® are the same as for Ihsternak foundation. Hov^ever, at z = 0, 
the shear force equation is to be modified. The boundary conditions, 
in this case, can be expressed as 


at 


as 


z = 0, = 0 j “ = IiaP 


ijL 

8z5 


,2 

A 0) 


n 


2EI 


Z-VOO y/oo,-|^^oo, eo , ^ 


82 


8Z 


(4.2.8) 


Crandell ( 14 ) derived an expression for the response of an infinite 
beam subjected to impulse load at its centre, resting on Wirkler 
foundation with dampiag included. He gave the expression for the 
boundary conditions aa described by (4.2.8), 
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4,2,2 General Solution Using Laplace Transform. 


The general solution of the homogeneous equation (4,2. l) based on 
equation (4.2,2) can be written as 

Kz) = + Cg ?2(z) + Cj ^^(z) + \{z) (4.2.9) 

where C , G , C , G are initial parameters, and $.(z) expression 
12^4 3 . 

are as follows 


4^(2) 


— — - (C^ co^ D z + cos C z) 
-2 ^ a a a a. ' 


-2 

C + D 
a a 


-2 -2 

C K ^ 

(z) = p ' ' sinh D z + sin C z) 

- + r D ® C ^ 

a a a a 


4_(z) = — s ; 3 p (cosh z - cos C z) 

^ G 4- D*^ • 

a a 


(4.2. 10) 


4^(2) = 


°a 


J (-L 

+ D 
a a 


sirih 


D^z - 


sin 


whei« the esxpzession for C and D are given by 

S . £L 

0^ = 4^ 

a 2 

^2 _ +/+/a^+ ^ 
a 2 

4 2 

ani = 4 (1 + 2DP + p ) 

Thus G„ and D are the functioneof ‘p*. 
a a ^ 


( 4 . 2 . 11 ) 


( 4 . 2 . 12 ) 
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(a) General Solutiori for Unite Beam oa lasternak PouMation. 

3Ihe general solution of equations (4.2,9) and (4,2, 10) can be solved 
for boundaiy conditions as described by equations (4.2.4), (4.2,5), 2he 
solution in la place transform for these boundary conditions can be 
written as 


y(z) =0^ fjCz) ~ ^4(2) . (4.2.13) 

yfoere C^, eire initial parameters and can be expressed as 

♦3'"(E) -4lJ(*j''(E) ♦^'(E) (E) 

ndi ) ? 3 ( e )»] 


a a' 


^ [i" (r) i^^’'’(R)- ?/'(R) ?3''’(R)- 4R^{|’3"(R) f/(R)- ?^"(R) ?3'(R) 

2 - -1 

l2'- -2Vl72 - <^ 1 "^®) <1>3(e)}} j 

^^a"" ^a^ ~ 


and 


[?4’'(R) f7'’(R)- ?^"(R) ?4''’(R)- 4 dJ{?4”(R) ^7’(r)- ?/'(R) $ 4 ’ (R) 




?,(«)- i,''(E) 


a a' 


?3''(e) i^"'(E)- f^"(E) «3 ’"(e)- A^l {Jj-Ce) 4/'(R)- t/lSlj'CE) 




a a' 


(4,2.14) 


Ihe solution at z = 0 can be written, after solving equations 

(4.2.13) and (4,2,14), as 
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D 


2(5 ^'^^(l+2Dp)(D cosh D E sin 0 R-C^siria D RccsCR) 

SL SL Q, Si 3* 3t 8,8. U.S. 


y(o,p) = Imp 


X^cd 


n 


-C D {2D^C^+ (^+^)cosh D E cos C„R 
aa aa^aa-^ a a 

+C D (D^-C^) sirii Ii„R sin 0„E} 
a a a a"^ a a J. 


2EI 


rc^E^(c^+D^){^cosh E^R sin OR + sirti D^R cos aR} 
L a a' a ^ a a a a a a 

-2W(i)l-ol)^'^^(cl+ D^)^(l+2Dp) cosh \B cos cj~j 

(4.2. 15) 

(b) Infinite Beam on Pasternak Foundation . As the response functions 

!D z 

at infinity have to be finite, the terms e a will have to be omitted in 
the expocession (4,2.15). After simplifying the expression, the general 
solution for deflection at the centre of the beam can be written as 


y(0,p) = Imp 


X uj 


n 


2EI 


r -5 R -D R ■ 

C e cos aR + E„e sin C„R 
a a a a, 

, 20 D (C^ + 5^) 

I- a a^ a a' 


(4.2.16) 


(c) Elnite Beam on Winkler Eoundation. If the shear coefficient A 
is made zero, the equation (4.2.15) will represent the solution of a 
finite beam on Winkler foundation, Eor boundaiy conditions described by 
(4.2.7), the general solution can be expressed as 


y(0,p) = Imp 


•i2 

X 


2EI 


1 + cos e^R cosh B^E 


sin B^R codi B^E + cos B^R sirii 3^® J 


(4.2.17) 


(d) Infinite Beam on Winkler foundation . For infinite beam, solution 
boils down to a simple expression, for the boundary conditions described 
by equation (4.2.8) and can be written as 
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2 

- ^ “n 1 

y(0,p) = Imp 3 

and general expression can be written as 


(4.2.18) 


,2 -e.z 
^ “n e ^ 

y(z) =r Imp^gj- — ^ — Qjos 0^z + sin g^z^ 


(4.2.19) 


4.3 imrERSIOIT OP LAEDACE TRAISPOM 


iaplace transform methods biingout a close correspondance between 
Tisco-elastic and elastic problems. In fact, ar^r problem of wave 
propagation in a visco-elastic medi-um can be reduced to Laplace inversion 
problem, so long as a solution can be found for an elastic solid subject 
to the same boundary conditions ( 7 ), Thus the techniques for inversion 
of LapLace transform can be applied for the present problem. Cost ( 13 ) 
has reviewed a number of inversion techniqpj.es. For impulse function, 

Post and Widder* s approximate solution ( 4 ) can be used for simple general 
equations. Other approximate solutions which can be used are Ter-haar's 
solution, Sdiapery's direqft method, Ihpoulis inversion method etc. (l5). 

In this section, solution for an infinite beam on Winkler foundation is 
derived using Itost & V/idder's approximate solution. As this method involves 
•Uie nth derivatives, it becomes quite tedious to develop expressions for 
finite beam on Winkler foundation or Pasternak foundation. Another 
inversion technique vhich m^ he used for complex problems, has been tried 
successfully for the inversion of Laplace transform expressions. This 
technique, developed by Berger (6) , is an inversion fomula for the 
LapLace transform in terms of a series expansion of Jacobi polynomials. 
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4,3.1. Jost & V/idder's Approximate Solution ( 4 ). 

She It>st and Widder' s approximate solution is used for the 
infinite beam on Winkler foundatioji. The expression is 


f(t) = limit 


n 


n+l jU 

(. 1 )“ i_f(p) 

nl , n 
dp 


p = r/t 


( 4 , 3 % 1 ) 


Ehe general solution for infinite beam using (4.3. 1 ) can be 
solved by breaking the trigonometric and exponential functions into 
series form, The general solution given by eqiiation (4.2,19) can 
be written in series form as 


(0 


y = Imp 


n 


2EI 


"1 2.1 .2^2 3 2.2 4 2.4 6 3 

Bl 1 


7ft4 2.8 8 5 2.16 _10 .7 , 

2 p ■i'-Trr z P, rrr z p, 


71 ~ ' 8*. “ "1 101 

facing the inversion term by term 


1 




-1 


1 


[574 


L 




-De it' A 


4 ,= I -’ r 4 - n = e— i _ (- pi .) 


rw 


e 




ire) 


For 0=5 


Ihe derivatives of the terms of the series are 


( 4 . 3 . 2 ) 


( 4 . 3 . 3 ) 


( 4 . 3 . 4 ) 


for 


= fT5/e + D)^ + (l-D^^ 


Y" =[(5/e tD)" t 


( 4 . 3 . 5 ) 



154 


9^3. 


9e' 


- — "'* Hg) ( N^)^~ ^( I^) 


(5) .5. 3. 5. 2.2/.. ^-714 


(]Sr^)“'^^'^(i5^)(l^) + iiL!|ilili2i2 (Np“'^^'^(ll2)(l^) 


- 5.5.3. 2. 2 (lip”^^^(llg) 


(4.3.6) 


v±iere = (|- + d)^ + s^ 

I 

2 6^ 




10 4. £- 

05 e4 


=4-^^ 

4 04 e5 


5 0^ ©5 




and s = 1-D 


Por 6 


5 = 1^/'* 
1 1 


5 3 


1 = . J M.isf.1.2 ( )-17/4( ,5^ 9,15]^ (N,)-«A(„ )3(. ) 

96 ^ ( 4 )^ ^ ^ ( 4 )^ ^ ^ ^ 


+ (5)^(3)^(U^)“^'^'^(N^)(l%) 
+ (5)^(3)^(N^)"‘^'^''’(l^)(Er2)-(5)^4.(3)^(Np”''/'^(Ng) 


.5/.s2 


(4.3.7) 
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itor 


8^ 

*^1 11. 


86' 


7»(|) , ^ , g A2 (u )~‘>5A(1L)5« 7 , . , (5) .5 , ;„ 2 x-1lA(jj.j3^^^j 

(4)^ ^ ^ (4)^ ^ ^ 

+ M!I3^(jjj- 7/4( ) (jj )2 

(4)^ ^ ^ ^ 




- (5)^3 

-(5)^.3. (2)^ (U^)'’/^(Kg) 


5 4 
8 3 


For 3^ = (N^) , 1 = -10.60.2 (Ug) 

80 


(4.3.8) 

(4.3.9) 


Denoting the derivatives as 

8^3 


A = -.l/5(5/0) -— , A = -l/5(5/0) — ^ ,4= -l/5(5/0) — ^ and 

2 ftA^ ^ ae^ ^ 30-^ 


88 ' 

5 5 

6 3 

Ag= -l/5(5/e) -'5^ 

86 


a5 r3 
.6 ^ ^ 

86 


6 9 


lEhe general solution can now be written as 


2 

/■ * N ^ *‘*n !” 

yC^js) = |^■|■■■ 


2 ^ 4 \ 6 1 , _7, ^ 8 - 

■ 2 -" ” 3 ^ 90 315 


■A^.-- Z-+ ^ Z-- ^ 4- ^ z ^ 


(4.3.10) 


4.3.2 lerlaee Inversion by Series Solution. 

ih,poulis ( 37 ) constructed an inversion formula of the one 
dimecBional lapiace transform in terms of series expansions of trigono- 
metric functions and Legendre Polynomials, Berger (6) constructed this 
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inversion, formula in teimis of series expansion of Jacobi polynond^s, 
Berger’s (6) technique seems to be better than the methods adopting 
Legendre polynomials or trigonometric functions. Berger’s technique 
is presented below in brief. 

The one dimensional Laplace transform of f ( 0 ) can be defined 


as 

00 

Kp) = / e”^®f(e) d0 

( 4 . 3 . 11 ) 

0 

Let f(x) = ^ a) (■^)^ 

( 4 . 3 . 12 ) 

p = (m+l)a + g-1 

( 4 . 3 . 13 ) 


vtoere a > 0, m is a positive integer and g is a parameter and 
assuming 




n^ 


the solution can be expressed as 


(4.3.14) 


P. P, . ^ (l-g/o) 

f(x) = i|j(x) (l-x) ^(1-Hx) ^ (-^) (4.3.15) 

(p.jPp)^ 

where is the Jacobi polynomial. 


Considering the time variable given by 

8 =[-i/^log (^) (4.3.16) 

and substituting (4.3.12), (4.3.13), (4.3.14), (4.3.15), (4.3.16) into 

( 4 . 5 . 11 ) ard changirg the order of summation, integration of equation 

( 4 . 3 . 11 ) gives 
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+1 ni+P- 


F [(231+1)0 + g -l”l=(l + a2™‘^^) I S J (l~x) ^(l-bc) 

L J nsO -1 

(Pl,P2)^ 




Squation (4.3.17) may be expressed as 

1 "2 


P (ni+l)c + g-1 


( p. + Po) m 


I g,E(a5n,p^,P2) 




l(p +m+l) I (p.+r.+l) for-m 

viiere B(m,n,p ,p ) = 7- ' " :i: 

nl |(-m) |Ip^P2+m+a+2) 

P*” denotes ganma f motion 


(4,3.17) 


(4.3.18) 


(4.3.19) 


Die coefficient may be computed from the equation (4.3.18) 


recursively and may be obtained as 

—I 

I 

k=0 




2a 


n 


(p.,>Pp) 


I (- 1 )^ C(k,n,p^,P2) P j^h+l)o + g-1^ (4.3.20) 


( P ^} P2 ) l ( p^+^-+1 frP2 '^“'^‘^ 


vshere h 


11 


(4.3.21) 


i(^) l(p,+P2+n+l) 

«j n 

and c(k,a,p^,P2) =^ (15.) (n+p^+Pg+l)-*- (n+P^+P2'^- (p^-^+l) (4.3.22) 


lEhe Jacobi Polynomial is given by 


n 


n 


(x) = I C(k,n,p ,p ) (^)^‘’ 
k=0 I a 


(4.3.23) 


Knoi’dng the values of Jacobi Polynomials for ary n, P^(p^jP2)» ard 
coefficient g^, the function 3|>(x) can be determined. !Baen the 
response function f(x) can be evaluated throu^ equation (4.3.15). 
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4.4 


RESDiTS 


Before deteimning the response values, the laplace inversion 
techuL^e was verified for displacement-tiine response at the centre 
of an infinite beam, resting on Winkler foundation. The exact solution 
in terms of fractional Bessel’s function, at the centre of an infinite 
beam, resting on ?/ihfcler four^ation, is available. The solution in 
laplace transform is given by (l4) 

Xu. . 


y(0,p) = ImPglj- 7— — rpV 

(1+2BP+P ) 


^2 
X u 


= Imp 


8® (p84-s2)5A 

2 2 2 2 
v^here s = -j-D P = (p+f) 


( 4 . 4 . 1 ) 


Thus after inversion 


y(0,e) = Imp -—y ( 7 — 


2;;:? 


1/4^ 


where 


W 

Imp =s — V 
^ g 0 


W v/~ - (14) 


Y = velocity of freely fallir^g mass © , just before 

O K 


( 4 . 4 . 2 ) 

( 4 . 4 . 3 ) 


striking the beaia. 


For the verification of the numerical inversion of Laplace 

transform, an impulse load (w = 40kg, v^* 441 cuv^sec.) was applied 

on an inanite beam of iiiss per unit length, 0,001 kg.sec.^cm^ and 

4 

momerd: of inertia 115.0 cm , resting on Winkler foundation 
2 

(K = 6000 kg/cm ). (Table 4. 1 ). For Laplace transform inversion the values 
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of constants for commutation are taken as = 3.0, 3*0 , 

0 = 50*0, g = 0.0. In figure 4.1 displaceroent time response is 
plotted by the dotted line for an infinite beam, resting on Wirfcler 
founda'cion, txsing numerical inversion technique. Solid curve shows 
the displacement time curve for beam, using the exact solution, 

liom this figure (lig. 4.l) it ca.a be observed that the response 
by numerical inversion is quite close to the response, by the exact 
solution. As the time is increasing, the time period of response, 
by numerical inversion, is increasing sli^tly. Also, the amplitudes 
are little hi^er compared to exact solution. For most of the 
practical problems, the response is important for initial two or three 
time cycles. In this range, the numerical resiflts can be seen to be 
matching closely with the exact solution (36). 

5ig, 4.2, shcnvs the displacement time response of an infinite 

beam on Winkler foundation for another set of soil and beam parameters 

(fable 4. l). In this case the spring constant of the soil medium is 
2 

800 ]s.g/om , She infinite beam of mass per unit length, equal to 
2 

.00016 — ^ , has the mcment of inertia 6.66 cm'". The impulse load 

was applied through a body of 10 kg, falli ng freely from the hei^t of 
100 cm, (w = 10kg, 441 cir/sec.). In this case also, the response 

can be seen to be close to the response, by ejact solution. With 
increasing tine, the response period is sli^tly increasing as was observed 
for the previous set of pjarameters, !ihe amplitude is also increasing 
but very marginally. For the initial two or three time cycles, the 
response is matching as in Fig. 4.1. 



EiS* 4.5 presents the peaks of the response curve with time 
for esact and numerical solution of displacement response at the 
centre of an infinite beam resting on Winkler foundation. VShile the 
peaks are occuilng at constant interval for the response by exact 
solution, the time period is increasing v/ith subsequent peaks, in 
response curve obtained by numerical inversion. Still, for initial 
two or three time cycles, it can be observed from the ligure (ilg. 4,5) 
that the time period is more or less constant. 

4.4,1 Infinite Beam 

In ELg, 4.4 the maximum displacement with time cycle is plotted 
tor various values of damping factor D, The solid lines daow the 
curves obtained by exact solution and dotted lines, by numerical 
inversion. It can be observed that with increase in time cycles, the 
maximum displacement is decreasing. Also, as the damping factor is 
increasing, the maxinmi displacement is decreasing. The rate of 
decrease in maximum displacement is hi^er for the damp»ed case, compared 
to undamped case. The maximum displacement values obtained by numerical 
inversion follows the same trend but the magnitude of displacement 
is little hi^er with increasing time, 

Elg, 4,5 presents the variation of maximum displacement with 
time for the ahother set of soil and beam parameters. The response 
curve is plotted frcra the results obtained by numerical inversion. As 
was observed in ilg, 4,4, the maximum displacement decreases with time. 

At lower time values, the rate of decrease is hi^er. ffith time 
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incregieing, the rate of decrease becomes less as in Hg. 4.4. Ihe 
maximum disyLacement values are laver for damped case and is decreasing 
with increase in damping factor, Ihe rate of decrease in maxinami 
displacement for damped case, is hi^er than, for the undamped case, 

Ihe variation of maxamum displacement of the first peak, with 
damping factor is shown in Fig, 4.6, for both the sets of soil and 
beam parameters, Sfexlmum displacement is decreasing mth increasing 
damping factor. 

Fig. 4.7 presents the variation of maximum displacemenb , along 
the length of the beam for the set 1 of soil and beam parameters given 
in table 4,1, Responses are plotted for three values of damping factor 
( 0 , 0 , 0 . 05 , 0 . 10 ), The maximum displacement (first peak) decreases 
along the length of the beam, ?liile at lower values of z, the rate of 
decrease is slow, the displacement decreases fast at hi^er values of z. 

For all the values of damping factors, the decrease in maximum displacement 
follows the same trend. 

The variation of maximum displacement along the length of the 
infinite beam is shown in. Fig. 4,8 for set 2 of soil and beam parameters 
(Table 4. I). Ihe maximum displacement is decreasing as 'z’ is increasing. 
For lov/er values of z, the rate of decrease is less. After that, the 
maximum displacement decreases fast. For different values of damping 
factor, the rate of decrease is nearly the same, as observed for set 1. 

Using Post & Widder's ( 4 ) approximate solution, the response 
curves for an infinite beam resting on Witifcler fourdation and subjected 
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to impulse load are drawn in ELg. 4.9 and 4.10. The curves are drawn 
for both damped and undamped cases. 

Tig. 4.9 shov/s the maodLimm dispLacement (displacement at the 
peak of first cycle) rarLation along the length of the beam for 
vLndamped and damped (D = 0.05) Winkler foundation, Ottie soil and beam 
parameters are as follows 

2 

Spring constant =5.0 kg/cm 

4 

Moment of inertia = 65 cm 

2 

Ivlass of the beam = .0016 

cm 

Impulse load = 40 kg sec. 

The maximum displacement is decreasing with the increase in z value. 

The maximxmi displacement is hi^er for ItO.O in comperison to the value 
of maximum displacement for D = 0.05. The rate of decrease for 
undamped and damped cases follow the same trend, i.e, at lower values 
of z the rate of decrease is slow, and at higher values of z, the 
decrease in maximum displacement is fast, 

Tig. 4.10 presents the variation of maximum displacement at the 
centre of the beam with damping factor (vdth the same data as given above). 
The maximum displac^aent is decreasing with the increase in damping 
factor. It can be observed that the decrease is aot quite linear. 

OSais curve is comiared v/ith the exact solution (l4). 

4.4.2 linlte Beam. 

For finite beam, resting on Winkler foundation, the response 
curves are drawn using Berger^ s (6) numerical inversion technique 
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Jig* 4.11 presents the variation of maximvuii displacement at the centre 
of the beam, vdLth the time, 5ig. 4,12 shoivs the variation of maxiinum 
displacement; along the length of the beam for damped and undamped cases. 

It can be observed from lig. 4.11 that maximum displacement: 
decreases fast as the time passes. As for infinite beam, the decrease 
in maximum displacement for damped foundation is rapid, that for 
undamped foundation. Again, the magnitude of maximum displacement of 
the finite beam at ary time t, is less than the amplitude for infinite 
beam at that time. In lig. 4.11, the variation of maximum displacement 
of finite beam is compiared with the displacement of infinite beam. 

Slg, 4,12, shows the variation of maximxmi displacement along the 
length of beam for damped and undamped cases. The maximum displacement 
along the length of beam for undamped foundation is higner than for 
I>=0,05. Ho?/ever, the rate of variation of maximum displacement along 
the length is more or less same for undamped and damped cases. 

Erom these results it can be observed that Beger's technique can 
be used successfully for the general solution of beams on elastic foun- 
dations subjected to impulse load. 
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Table 4.1 

Soil and Beam Parameters used for Ilmerical Results by Inversion 

Technique 


Set 

Spring 

constant 

kg/cm^ 

Modulus 
of Ela- 
sticity 
kg^cm 

Moment 

of 

Inertia 

4 

cm 

Ivlass of 

the Beam 

, 2 
kg. sec 

2 

cm 

VVeighf of 
falling 
body kg. 

Hei^t of 
Ehll cm. 

Set Ko.l 

6000.00 

2100000.0 

113.0 

0.001 

40.00 

100.00 


Set I0.2 800 .00 2100000.0 


6.66 0.00016 


10.00 


100.00 
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Fig. 4.1 Variation of displacement with time of an infinite 
beam on Winkler foundation subjected to impulse 
load (set 1) 



Fig. 4.2 Variation of displacement with time of an infinite 


beam on Winkler foundation subjected to impulse 
load (set 2) 
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Fig. 4.4 Variation of maximum displacement Fig.4.3 Comparison of time period of 

at the centre of the beam with , the response by numerical 

time (set 1). inversion with that of exact 

solution. 
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Fig. 4-6 Variation of maximum displacement Fig.4.5 Variation of maximum displacement at 
with damping factor (Bergers the centre of the beam with time 

techhique) (Bergers technique) (set 2) 


Maximum displacement (in mm) 



0.0 0.1 


0.2 0.3 0.4 0.5 


Distance from origin (x/l ) 


Fig. 4.7 Variation of maximum displacement along 
the length of beam (Berger's technique) 
(set 1) 
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Fig. 4.8 Variation of maximum displacement along 
the length of beam (Berger's technique) 
(set 2) 
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Fig. 4.9 Variation of maximum displacement 
along the length of beam. 

(Post Ut V\/idder's technique) 
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Fig. 4.10 Variation of maximum displacement at 
the centre of beam with damping 
factor (Post & Widders technique) 
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at the centre of the finite beam of finite beam (R=0.6) along the 

^(RrO.6) (Bergers technique) length of the beam (Berger's 

technique) 



CHAPTER 5 


CONCLUSIORS MD RECOIlEKD^iTIONS 

5.1 CORCLUSIOtlS 

Ex-DeriTiBntal and analytical investigations have been carried out 
on some problems in the area of dynamic soil-struoture_interactionsin 
the present work. Results have been discussed and conclusions have 
been presented in each chapter separately. However, general conclusions 
are presented belov/ from the present investigations. 

‘ To facilitate the correct evaluation of soil parameters which 
characterize the behaviour of soil medium, studies have been made on 
the effect of magnitude of force levels and flexibility of the footing 
on the foundation response . Rrom the experimental results it can 
be observed that with the increase in force level, the magnitude of 
displacement increases more or less proportionately. The theoretical 
analysis of the footings on soil medium, shows that the displacement 
of the footing is directly proportional to the force applied. The 
flexibility of the footing has very little effect on spring constant. 
iSnalytical solutions based on elastic theory show that there is no 
difference in spring constant values due to variation in flexibility. 

The studies on the effect of size and shape of the foundation on soil 
parameters reveal that with length and width variation, the soil 
parameters such as spring constant, damping factor etc. vary considerably. 
While the spring constant is increasing with increasing l/b ratio, the 
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damping factor is decreasing. Thus, in the analysis of rectangular plates 
of higher l/b ratios and long heams, these soil parameters may have to 
he modified. 

From the anal5rtical results of the beam on elastic foundations 
(Winkler model and Pasternak model) carrying a concentrated mass on it 
and subjected to steady state dynamic force,it oanbe seen that the effect 
of concentrated mass on the displacement response is quite pronounced. 

The maximum displacement decreases very fast as the concentrated mass 
increases on the beam. However, for the beams resting on Pasternak 
foundation, the shear coefficient of the shear layer of the foundation 
model, plays a more important role in decreasing the displacement, not 
only at the centre, but also, along the length of the beam. For certain 
values of shear coefficient, the effect of concentrated mass is not 
noticeable. To compare the theoretical results, experiments have been 
conducted on a long beam carrying a concentrated mass and resting on 
the natural soil deposit . These experimental values are in close agreement 
with the theoretical results. 

The beam resting on elastic foundation has again, been analysed 
for the transient load (impulse load). The beam subjected to impulse 
load has been analysed using Laplace transform. Tne solutions have been 
obtained using numerical inversion of Laplace transform. While the 
Post and Widder’s formula for numerical inversion was used for infinite 
beam on elastic foundation, the Berger's Laplace inversion technique 
was used for finite and infinite beams on elastic foundations. However, 
Berger's technique has been shown to be wore useful for complex problems. 
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while Post and Widder's formula can he used for sin^ile cases. These 
numerical results show that the displacement along the length of the 
beam decreases as the distance from origin increases. With damping of 
the foundation, also, the displacement decreases. 

5,2 HECOI.tEHDATIONS 

The following topics are recommended for further stuies 

(1) Experimental investigations on the effect of geometry on prototype 
foundations. 

(2) Experimental investigations for contact pressure distributions 
for dynamic loads. 

(3) Effect of shear deformation and rotary inertia of the beam on the 
beam-foundation response. 

(4) Beams subjected to periodic and aperiodic loads such as triangular 
pulse, rectangular pulse etc. 
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